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CHAPTER  I 


INTRODUCTION 


Several  Radar  Target  Identification  (RTI)  techniques  have  been  developed  at 
The  Ohio  State  University  in  recent  years.  Using  the  ElectroScience  Laboratory 
compact  range  a  large  data  base  of  coherent  RCS  measurement  has  been  con¬ 
structed  for  several  types  of  targets  (aircraft,  ships,  and  ground  vehicles)  at  a 
variety  of  polarizations,  aspect  angles,  and  frequency  bands.  This  extensive  data 
base  has  been  used  to  analyze  the  performance  of  several  different  classification 
algorithms  through  the  use  of  computer  simulations. 

In  order  to  optimize  classification  performance,  Ksienski  [l]  concluded  the 

/' 

radar  frequency  range  should  lie  in  the  Rayleigh-resonance  frequency  range,  where 
the  wavelength  is  on  the  order  of  or  larger  than  the  target  size.  For  aircraft 
and  ships  with  general  dimensions  on  the  order  of  10  meters  to  100  meters  it 
is  apparent  that  the  High  Frequency  (HF)  band  provides  optimal  classification 
performance.  Since  existing  HF  radars  are  currently  being  used  for  detection  and 
tracking  of  aircraft  and  ships  of  these  dimensions,  it  is  natural  to  further  investigate 
the  possibility  of  using  these  existing  radars  as  the  measurement  devices  in  a  radar 
target  classification  system.  — —  (  ,  N  *  r'  '  ''f 

A  general  radar  target  classification  system  is  shown  in  Figure  1.  The  external 
environment  consists  of  a  target  which  is  to  be  classified,  clutter  source(s),  noise 
source(s),  and  possible  calibration  reference(s).  The  environment  is  observed  by 
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the  inultifrequency  radar  through  the  propagation  media.  The  signal  processor 
provides  radar  data  which,  when  combined  with  a  priori  knowledge  of  the  external 
environment  and  knowledge  of  the  radar  system  and  propagation  media,  allows  the 
estimation  of  target  features  to  be  used  for  classification.  These  estimated  features 
of  the  unknown  target  are  compared  with  high  accuracy  features  of  several  targets 
in  a  catalog  set.  The  classification  algorithm  provides  a  measure  of  the  similarity 
of  the  unknown  target  to  each  catalog  target.  The  catalog  target  whose  features 
most  closely  represent  that  of  the  unknown  targets’  features  is  chosen  as  the  type 
of  the  unknown  target. 

This  study  investigated  various  propagation  conditions  which  allow  different 
target  features  to  be  estimated.  These  conditions  were  synthesized  by  various 
channel  models.  The  incorporation  of  these  channel  models  into  a  radar  target 
identification  system  computer  simulation  provided  estimates  of  resulting  classifi¬ 
cation  performance.  As  a  result,  the  relationship  between  channel  conditions  and 
classification  performance  were  found. 

IfF  radars  are  generally  catagorized  by  the  propagation  mode  by  which  they 
are  designed  to  operate,  either  1)  surface  wave  or  2)  skywave.  Each  propagation 
mode  incites  different  characteristics  into  the  classification  system,  and  therefore 
should  be  investigated  independently. 

Because  of  their  unique  propagation  mechanism,  surface  wave  radars  have  the 
ability  to  look  beyond  the  horizon  with  detection  ranges  up  to  300  km.  Surface 
wave  attenuation  is  a  function  of  many  parameters,  most  importantly  ground  con¬ 
duction  and  frequency  [6].  In  order  to  achieve  practical  detection  ranges  for  given 
transmitted  power  levels  the  radar  wave  usually  must  propagate  over  water.  Fur¬ 
thermore  only  a  vertically  polarized  wave  propagates  significantly,  thus  limiting 
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CHAPTER  I 
INTRODUCTION 

Several  Radar  Target  Identification  (RTI)  techniques  have  been  developed  at 
The  Ohio  State  University  in  recent  years.  Using  the  ElectroScience  Laboratory 
compact  range  a  large  data  base  of  coherent  RCS  measurement  has  been  con¬ 
structed  for  several  types  of  targets  (aircraft,  ships,  and  ground  vehicles)  at  a 
variety  of  polarizations,  aspect  angles,  and  frequency  bands.  This  extensive  data 
base  has  been  used  to  analyze  the  performance  of  several  different  classification 
algorithms  through  the  use  of  computer  simulations. 

In  order  to  optimize  classification  performance,  Ksienski  [1]  concluded  the 
radar  frequency  range  should  lie  in  the  Rayleigh-resonance  frequency  range,  where 
the  wavelength  is  on  the  order  of  or  larger  than  the  target  size.  For  aircraft 
and  ships  with  general  dimensions  on  the  order  of  10  meters  to  100  meters  it 
is  apparent  that  the  High  Frequency  (HF)  band  provides  optimal  classification 
performance.  Since  existing  HF  radars  are  currently  being  used  for  detection  and 
tracking  of  aircraft  and  ships  of  these  dimensions,  it  is  nat  ural  to  furt  her  investigate 
the  possibility  of  using  these  existing  radars  as  the  measurement  devices  in  a  radar 
target  classification  system. 

A  general  radar  target  classification  system  is  shown  in  Figure  1.  The  external 
environment  consists  of  a  target  which  is  to  be  classified,  clutter  sotirce(s),  noise 
source(s),  and  possible  calibration  reference(s).  The  environment  is  observed  by 
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the  multifrequency  radar  through  the  propagation  media.  The  signal  processor 
provides  radar  data  which,  when  combined  with  a  priori  knowledge  of  the  external 
environment  and  knowledge  of  the  radar  system  and  propagation  media,  allows  the 
estimation  of  target  features  to  be  used  for  classification.  These  estimated  features 
of  the  unknown  target  are  compared  with  high  accuracy  features  of  several  targets 
in  a  catalog  set.  The  classification  algorithm  provides  a  measure  of  the  similarity 
of  the  unknown  target  to  each  catalog  target.  The  catalog  target  whose  features 
most  closely  represent  that  of  the  unknown  targets’  features  is  chosen  as  the  type 
of  the  unknown  target. 

This  study  investigated  various  propagation  conditions  which  allow  different 
target  features  to  be  estimated.  These  conditions  were  synthesized  by  various 
channel  models.  The  incorporation  of  these  channel  models  into  a  radar  target 
identification  system  computer  simulation  provided  estimates  of  resulting  classifi¬ 
cation  performance.  As  a  result,  the  relationship  between  channel  conditions  and 
classification  performance  were  found. 

HF  radars  are  generally  catagorized  by  the  propagation  mode  by  which  they 
are  designed  to  operate,  either  1)  surface  wave  or  2)  sky  wave.  Each  propagation 
mode  incites  different  characteristics  into  the  classification  system,  and  therefore 
should  be  investigated  independently. 

Because  of  their  unique  propagation  mechanism,  surface  wave  radars  have  the 
ability  to  look  beyond  the  horizon  with  detection  ranges  up  to  300  km.  Surface 
wave  attenuation  is  a  function  of  many  parameters,  most  importantly  ground  con¬ 
duction  and  frequency  [6].  In  order  to  achieve  practical  detection  ranges  for  given 
transmitted  power  levels  the  radar  wave  usually  must  propagate  over  water.  Fur¬ 
thermore  only  a  vertically  polarized  wave  propagates  significantly,  thus  limiting 
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Figure  1:  General  classification  system. 


the  polarization  diversity  of  our  measurement  system. 

HF  skywave  radar  utilizes  ionospheric  reflection  as  the  propagation  mecha¬ 
nism.  The  skywave  radar  has  an  average  detection  range  of  1000  km  to  4000  km, 
which  when  coupled  with  an  antenna  steerable  in  azimuth  up  to  120°,  provides 
coverage  of  over  15  million  square  kilometers  [12].  The  ionosphere  supports  HF 
radio  propagation  at  altitudes  of  100  to  500  km  and  consists  of  ionized  regions 
referred  to  as  layers  which  are  commonly  labeled  D ,  E,  F\ ,  and  Fo-  Ions  at  these 
altitudes  are  mainly  produced  by  solar  radiation  which  includes 'particle  radiation, 
ultraviolet  light,  and  x-rays.  In  general,  it  is  the  electron  density  distribution  at  a 
given  altitude  which  effects  the  propagation  of  the  EM  wave  to  the  greatest  extent. 

The  actual  ion  distribution  at  any  time  is  a  function  of  many  factors  includ¬ 
ing  source  radiation  intensity,  angle  between  zenith  and  solar  radiation  rays,  ion 
distributions,  and  various  ion  combination  mechanisms  as  described  in  Davies  [14]. 
Therefore  the  ionosphere  is  continually  changing  and  must  be  monitored  in  real¬ 
time  to  accurately  determine  its  present  structure.  It  is  the  variation  in  electron 
density  with  altitude  which  causes  variation  in  the  refractive  index  and  results 
in  the  bending  of  the  ray  path  and  its  return  to  earth.  It  is  important  to  note 
the  ionosphere  only  supports  propagation  over  a  limited  band  of  frequencies,  and 
signals  of  different,  frequencies  may  travel  along  different  ray  paths. 

The  earth’s  magnetic  field  also  affects  the  propagation  of  the  radar  wave.  The 
magnetic  field  and  ionosphere  produce  a  magnetoionic  medium  which  supports 
propagation  modes  with  specific  polarizations.  These  polarizations  are  a  function 
of  the  magnetic  field,  parameters  of  the  ionosphere,  and  the  direction  of  propaga¬ 
tion  as  given  in  Davies  [14] .  As  a  result,  a  transmitted  wave  of  a  given  polarization 
may  split  into  two  ray  paths  which  may  or  may  not  travel  the  same  route  to 
the  target.  Therefore,  in  general,  the  polarization  of  the  transmitted  wave  is  not 


necessarily  the  polarization  of  the  wave  which  illuminates  the  target,  and  the  po¬ 
larization  of  the  field  scattered  from  the  target  is  not  necessarily  the  polarization 
of  the  return  wave  at  the  receive  antenna.  Researchers  such  as  Pilon  and  Headrick 
[24]  have  addressed  this  situation  as  it  applies  to  the  RC'S  estimation  problem. 
This  phenomena  must  be  carefully  investigated  when  selecting  the  polarization  of 
the  catalog  sets. 
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CHAPTER  II 


HF  RADAR  SIGNAL  PROCESSING 


2.1  DOPPLER  PROCESSING 


Many  operation  HF  radars  are  used  for  the  detection  and  tracking  of  targets 


on  or  over  ocean  waters.  A  typical  radar  cell  illuminated  by  a  surface  wave  radar  is 


shown  in  Figure  2.  The  width  of  the  observation  cell  is  determined  by  the  antenna 


beamwidth  at  the  radar  operating  frequency,  while  the  length  is  determined  by 


the  time  gating  of  the  received  signal.  To  observe  a  cell  at  a  distance  R\,  from  the 


antenna  with  a  cell  length  A R,  the  received  signal  is  gated  in  time  with  values: 


r2  =  rj  + 


where  c  is  the  speed  of  light,  and  (r2  —  t\)  is  the  time  gate  used.  With  a  multi¬ 


frequency  radar,  the  beamwidth  of  the  antenna  may  change  since  the  radiation 


pattern  is  a  function  of  frequency.  Therefore  the  size  of  the  cell  varies  with  changes 


of  the  radar’s  operating  frequency. 


In  operation,  the  radar  transmits  a  continuum  of  modulated  pulses  at  some 


given  pulse  repitition  frequency  (PRF).  A  total  of  N  pulses  are  t  ransmitted  during 


the  coherent  integration  time  To  where: 


N  =  PRF  ■  Tc 


.<3  vf»^-£k  **  til 


> 


:  w  vy  wvwuit/wvusvtwvvvvvvvv^- vvir.nrairj  i 


Therefore  N  coherent  data  records  are  produced  from  the  energy  scattered 
from  the  radar  cell.  These  data  are  Doppler  processed  to  produce  a  power  spectral 
density  estimate  containing  N  data,  with  a  frequency  resolution  J-  and  a  band- 

c 

width  from  —  to  +  — ^ .  Theoretically,  greater  resolution  may  be  obtained 
by  increasing  of  the  coherent  integration  time;  however  changes  in  the  ionospheric 
ray  path,  sea  state,  and  target  place  a  limit  on  the  practical  integration  time  [22]. 
The  resulting  spectra  shown  in  Figure  3  exhibits  1 )  a  return  from  the  target  Pj{f ) 
with  induced  doppler  shift  /j,  2)  two  large  clutter  returns  P+g(f)  an^  P- #(/)  at 
\fg  and  —fg  respectively,  and  3)  other  relatively  low  power  clutter  covering  the 
spectrum. 

2.2  DOPPLER  SPECTRAL  ESTIMATE  INFORMATION  CONTENT 

Information  contained  in  the  spectral  estimate  may  be  used  to  estimate  a 
number  of  target  features  for  classification  purposes.  Using  fj,  the  relative  velocity 
of  the  target  may  easily  be  estimated,  giving  the  RTI  designer  knowledge  of  the 
probable  class  of  the  target  (ships,  general  aircraft,  high  speed  aircraft,  etc.).  In 
addition,  estimates  of  the  power  returned  from  the  target  Pj(f)  and  the  power 
returned  from  some  reference  Pg(f)  may  be  used  to  estimate  RCS  related  target 
features. 

The  two  large  power  densities  at  -\-fg  and  —  fg  result  from  electromagnetic 
wave  scattering  from  ocean  wave  sets  of  a  certain  wavelength.  The  ocean  waves 
which  excite  this  Bragg  scattering  mechanism  have  a  wavelength  dependent  upon 
the  radar  frequency  [15].  As  explained  in  detail  by  Trisna  [15],  Barrick  [5],  and 
others,  the  amplitude  of  the  Bragg  returns  are  largely  dependent  upon  sea  state. 

Large  variations  in  the  estimates  P+g(f)  «md  P-Bif)  about  some  means 
( P+g(f ))  and  {P~g(f))  may  be  observed  over  a  given  observation  period  as  noted 
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by  Headrick  and  Pilon  [24].  Therefore  the  Bragg  returns  may  best  be  described 


as  a  random  process  with  its  mean  being  a  function  of  the  sea  state  at  that  time. 


An  understanding  of  this  process  is  important  if  the  Bragg  lines  are  to  be  used  as 


a  reference  in  the  target  feature  estimation  process. 


2.3  STATISTICS  OF  SPECTRAL  COMPONENTS 


As  discussed,  the  height  of  a  target  return  or  Bragg  line  in  a  given  doppler 


spectra  is  only  an  estimate  of  the  power  returned  from  the  respective  scatter. 


Consider  creating  (Af  —  1)  additional  doppler  PSD  estimates  of  the  same  cell  to 


give  a  total  of  M  doppler  PSD  estimates.  Furthermore,  create  these  additional 


estimates  by  observing  the  cell  over  a  continuous  time  span  IqBS »  where 


tQBS  —  ^EST  x  M 


as  shown  in  Figure  4.  Examples  of  real  data  obtained  in  this  manner  may  be  found 


in  Pilon  and  Headrick  [24]. 


Note  in  Figure  4  that  each  spectra  gives  a  different,  estimate  for  each  Bragg 


peak  as  well  as  the  target  spectral  peak.  However,  by  averaging  M  samples  taken 


over  M  consecutive  CIT  periods  (k  =  1,...M)  one  may  produce  an  averaged  doppler 


PSD  as  shown  in  Figure  4.  The  averaged  estimated  power  returned  from  the  target, 


approaching  Bragg  wave  set,  and  receding  Bragg  wave  set  are  given  by  Eq.  (2.5), 


Eq.  (2.6),  and  Eq.  (2.7)  respectively. 


,  M 
1  ST  kb 


Wl)=77E*Wl) 


1  M 

P+B(fl)=  T7  Ekp+B(fl) 


i  M 

P-bUi)=  T7 


i Vv i.v Vf 


If  each  of  the  samples  on  the  right  hand  side  of  Eq.  (2.5)  are  independent  and 
identically  distributed,  the  application  of  the  central  limit  theorem  as  discussed  in 
Degroot  [10]  gives: 


lim 
M  — >oo 


i  M 

M  £  kprU\) 

k- 1 


Jim  =  (PT(fl)) 

M  — »oo 


(2.8) 


where  (Pj^if i))  is  the  expected  value  of  the  returned  power  from  the  target  as 
discussed  earlier.  Likewise,  application  on  the  central  limit  theorem  to  Eq.  (2.6) 
and  Eq.  (2.7)  yields  Eq.  (2.9)  and  Eq.  (2.10)  respectively. 


lim 

J1/-.00 


1  M 
M  k  =  1 


Jim  P?B(h)  =  (P+B(h)) 

M  — »oc 


(2.9) 


lim 

M  — >oo 


1 

M 


M 


k=  1 


Jim  P*!g(fi)  =  (P-B(fi)) 

M  — *oo 


(2.10) 


Although  this  implies  an  infinite  observation  time  which  is  impractical,  it  also 
implies  that  we  can  continually  improve  the  accuracy  of  our  estimates  by  taking 
more  samples. 

Both  the  CIT  and  the  observation  time  must  be  carefully  chosen  in  order  to 
satisfy  the  requirement  of  independence  and  identically  distributed  samples.  The 
observation  time  must  be  selected  to  assure  the  identical  distribution  of  samples.  A 
significant  change  in  target  flight  path  during  the  observation  time  may  cause  the 
observed  RCS  to  change  and  therefore  {Pj{f \))  changes  voiding  the  identically 
distributed  criteria.  Likewise  a  change  in  sea  state  during  the  observation  time 
may  alter  (P+g(f  j ))  and  {P~g(f  1 ))  with  the  same  result  .  In  general,  however,  sea 
states  may  remain  unchanged  for  hours  [15].  In  addition,  significant  changes  in  the 
observation  path  (such  as  passage  of  ionospheric  disturbance  for  sky  wave  radars) 
limits  the  observation  time.  In  fact,  fluctuations  in  the  ionosphere  generally  impose 
the  greatest  restriction  on  the  maximum  observation  time. 
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Figure  4:  Averaging  of  independent  spectral  estimates. 


There  is  a  limit  to  the  shortest  possible  integration  time  that  produces  inde¬ 
pendent  samples  of  the  Bragg  spectral  lines.  An  experimental  study  by  Barrick 
[13]  concluded  that  an  integration  time  of  25  seconds  or  greater  insures  indepen¬ 
dence  of  spectra,  irregardless  of  sea  state  and  cell  size,  at  frequencies  in  the  HF 
band  above  7  MHz. 

The  remaining  question  is,  given  a  CIT  of  25  seconds,  what  is  the  relationship 
between  the  standard  deviation  values  of  Pj?  (fl  )?  and  P+g{f\ ),  P^g(f\ ),  and  Ml 
As  discussed  by  Barrick  [13],  without  assuming  Gaussian  statistics,  there  is  little 
to  do  other  that  record  massive  amounts  of  data  and  extract  empirical  results. 
Using  a  Gaussian  assumptions  Barrick  derived: 

„  (<(/)) 

v^7 

where  a ( ,  s  the  standard  deviation  of  the  random  variables  given  in  Eq.  (2.6) 
*  B  w ) 

and  Eq.  (2.7).  Figure  5,  from  Barrick  [13],  shows  the  90%  confidence  interval  lines 
derived  using  the  same  statistical  assumption.  These  curves  give  the  RTI  designer 
insight  into  the  tradeoffs  between  the  observation  time  and  the  accuracy  of  the 
estimates  P+g(f\)  and  P^g(f\  )•  This  is  important  if  P^g(f\)  or  P^g(f\ )  are  to 
be  used  in  the  estimation  of  target  features  for  classification. 

2.4  FREQUENCY  SELECTIVE  FADING 

The  Doppler  spectral  estimate  resulting  from  data  accumulated  via  sky  wave 
propagation  is  distored  to  varying  degrees  by  the  ionospheric  path.  The  distortion, 
referred  to  as  ionospheric  contamination,  is  a  function  of  parameters  that  may  or 
may  not  be  under  the  control  of  the  radar  operator.  It  is  imparitive,  if  data  from 
the  resulting  spectra  are  used  in  an  estimation  and  classification  procedure,  that 
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Figure  5:  Relationship  between  deviation  of  average  Bragg  line  amplitude  and 
number  of  independent  samples  used  for  averaging. 
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(the  CIT).  The  resulting  spectra  will  be  offset  in  floppier  frequency  by: 
,  _  2vj  2vj 


(2.12) 


where  t'}  is  the  relative  ionospheric  velocity  with  respect  to  direction  of  radar  wave 
propagation  during  the  pass  and  A#  is  the  radar  wavelength.  This  phenomena, 
referred  to  as  Doppler  shifting,  is  easily  detected  in  the  skewing  of  the  spectra  as 
shown  in  Figure  6,  and  can  be  easily  compensated  for. 

Next,  consider  the  case  where  the  C-IT  is  relatively  long  such  that  the  relative 
velocity  of  the  propagating  layer  has  undergone  considerable  change.  Each  pass¬ 
ing  wave  may  have  a  slightly  different  doppler  component  induced  by  the  quickly 
changing  ionosphere.  The  resulting  contamination  of  the  spectra,  known  as  broad¬ 
ening,  is  shown  in  Figure  7. 

Spectral  broadening  is  highly  undesired  for  several  reasons.  First,  any  reduc¬ 
tion  in  the  amplitude  of  the  target  ret  urn  will  make  the  ret  urn  more  succeptable  to 
other  noise,  thus  making  any  target  feature  estimated  from  the  data  less  reliable. 
Second,  the  return  from  the  target  may  not  be  significantly  greater  than  the  higher 
order  sea  return  and  noise  surrounding  it,  and  thus  the  target  may  be  easily  lost 
in  the  clutter  surrounding  it.  Third,  if  the  target  has  a  Doppler  frequency  which  is 


1 


the  distortion  be  removed  or  at  least  understood  so  it  may  be  compensated  for  in 
the  process. 

Consider  a  target  being  illuminated  by  a  single  wave  which  travels  along  a 
single  ionospheric  path.  The  ionosphere,  which  is  in  constant  motion,  induces  a 
Doppler  shift  on  the  rate  and  direction  of  the  ionospheric  movement  .  Now  consider 
illuminating  this  target  N  times  by  waves  which  propagate  by  the  same  layer  whose 
relative  velocity  remained  unchanged  during  the  time  the  data  was  accumulated 


v-vyv 


w: 
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near  a  Bragg  frequency,  the  resultant  spreading  of  the  Bragg  line  may  essentially 
cover  the  target. 

An  additional  cause  of  Doppler  spreading  is  multipath  propagation  [22].  The 
path  or  paths  the  wave  takes  between  the  transmit  antenna  and  the  target  is  a 
function  of  many  radar  parameters  including  frequency,  launch  angle,  and  polar¬ 
ization,  as  well  as  real-time  ionospheric  conditions.  Different  ionoshperic  layers 
may  support  propagation  of  waves  of  the  same  frequency.  This  may  result  in  two 
or  more  waves  with  different  induced  doppler  shifts  both  containing  information 
about  the  target  being  processed  together,  with  a  resulting  spectra  which  exhibits 
spreading. 

A  third  cause  of  ionospheric  distortion  results  from  horizontal  spacial  inho- 
mogenity  within  a  given  layer  [22].  A  ray  path  utilizes  an  area  of  a  given  layer 
proportional  in  size  to  the  beamwidth  of  the  antenna.  The  larger  the  beamwidth, 
the  larger  the  area  used  for  refraction,  and  the  greater  the  effect  of  the  spreading  as 
investigated  by  George  and  Marseca  [22].  Since  beamwidth  scales  inversely  with 
frequency,  we  recognize  that  the  effect  of  spacial  inhomogenity  will  decrease  at 
higher  frequencies  for  propagation  via  a  given  layer.  The  RTI  system  designer  is 
interested  in  an  antenna  with  a  minimal  beamwidth;  however  physical  size  poses 
a  design  restriction.  For  example,  a  beamwidth  of  1°  may  require  antenna  dimen¬ 
sions  on  the  order  of  2  km  [12], 

In  order  to  use  doppler  spectral  data  in  an  estimation  process,  one  first  desires 
a  measure  of  the  degree  of  the  spectral  contamination.  In  general,  the  width  of  the 
Bragg  lines  gives  a  reasonable  measure  of  ionospheric  contamination  [22].  A  great 
deal  of  research  has  gone  into  this  area  as  well  as  the  development  of  processes 
designed  to  clean  up  the  distorted  spectra  [22]. 

An  alternate  approach  is  to  minimize  distortion  by  judicious  select  ion  of  radar 


parameters.  The  selection  of  an  observation  time  will  result  in  a  trade-off  between 
resolution  and  broadening.  Since  spectral  peaks  are  largely  uncorrelated  after  25 
seconds  (see  Section  2.3)  we  desire  the  shortest  CIT  above  25  seconds  that  meets 
the  minimum  spectral  resolution  demands  of  the  system.  This  allows  the  greatest 
number  of  independent  spectra  to  be  amassed  in  the  least  total  observation  time. 
Coherent,  integration  times  ranging  from  25  seconds  to  100  seconds  are  common. 
The  amount  of  disturbance  a  given  ionospheric  region  is  undergoing  determines  the 
degree  of  spreading  for  a  given  CIT.  Typically,  the  E  and  E$  regions  are  considered 
most  stable  but  do  not  always  exist.  During  sunrise  and  sunset,  the  ionospheric 
regions  undergo  rapid  change  and  the  degree  of  broadening  significant  ly  increases. 

The  judicious  choice  of  radar  operating  frequencies  can  also  greatly  reduce 
the  effect  of  ionospheric  contamination.  There  are  however,  important  trade-offs 
between  the  number  of  frequencies  over  which  one  may  obtain  data  and  the  quality 
of  that  data.  The  RTI  system  designer  desires  estimates  of  feat  ures  at  a  number 
of  frequencies.  Through  the  use  of  vertical  incidence  ionograms  and  sweep  of 
frequency  backscatter  echoes  one  can  determine  the  frequency  bands  least  likely 
to  result  in  multipath  distortion  [25]. 

One  possible  method  the  RTI  system  designer  may  choose  is  to  use  an  ex¬ 
isting  measure  of  spectral  contamination  with  a  limit,  throwing  away  all  spectra 
contaminated  beyond  that  limit,  and  making  his  estimates  and  classification  de¬ 
cision  based  only  on  the  remaining  data.  This  method  is  introduced  by  Georges 
in  [22].  A  second  approach  is  to  create  feature  estimates  using  all  possible  spectra 
and  then  weight  the  features  in  the  classificaton  process  according  to  the  degree 
of  distortion  of  the  data  from  which  the  features  were  taken.  A  third  approach 
is  to  estimate  features  from  the  doppler  spectra  whose  values  are  not  effected  by 
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the  spreading.  All  major  peaks  in  the  spectra  are  equally  effected  by  spreading 
such  that  the  ratios  ’  an<^  re,na’n  unchanged  [23] 

Therefore  use  of  these  ratio  in  feature  estimation  may  provide  estimates  which  are 
relatively  unaffected  by  Doppler  spreading. 
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CHAPTER  III 


MULTI-FREQUENCY  PARAMETER  ESTIMATION 


3.1  GENERAL  CONSIDERATIONS 


All  RTI  algorithms  utilize  estimates  of  various  features  of  the  unknown  target 
to  classify  the  target.  Features  such  as  radar  cross  section,  intrinsic  target  phase, 
and  target  speed  may  be  considered  good  target  descriptors.  However,  there  is  a 
trade-ofT  between  the  number  of  features  used  in  a  system,  the  cost  of  the  system 
(processing  time,  etc.),  and  the  resulting  classification  performance.  As  a  result, 
research  has  been  directed  at  finding  the  optimal  types  and  numbers  of  features 
which  result  in  optimal  performance  at  minimal  cost  given  little  or  no  constraints 
in  our  ability  to  obtain  estimates  of  these  features  [2],  [4]  and  [8].  In  an  operational 
HF  radar  classification  system  the  numbers  and  types  of  estimated  features  may 
be  severely  limited  due  to  the  constraints  imposed  by  the  HF  channel. 

The  purpose  of  this  section  is  to  study  the  constraints  the  HF  radar  system  im¬ 
poses  on  the  classification  system.  Given  these  constraints,  we  can  then  investigate 
which  and  how  many  of  the  features  available  should  be  used  to  obtain  optimal 
performance.  Previous  research  has  led  to  the  development  of  RTI  algorithms 
which  incorporate  multi -frequency  measurements  of  RCS  and  intrinsic  phase  in 
various  forms  as  classification  features.  The  number  and  selection  of  frequencies 
to  be  used  has  been  studied  with  regard  to  resulting  classification  performance  [2]. 
Different,  classification  algorithms  may  place  additional  restrictions  on  frequency 
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selection.  However,  in  an  operational  HF  RTI  system,  the  bandwidth  is  likely  to  be 
constrained  by  propagation  parameters,  thereby  constraining  the  number  and  type 
of  features  estimated.  The  different  propagation  modes  of  surface  wave  and  sky 
wave  radars  result  in  different  types  of  bandwidth  restraints  on  the  two  systems, 
therefore  they  are  investigated  independently. 

3.2  FREQUENCY  BAND  CONSIDERATIONS 

The  frequency  range  of  a  surface  wave  radar  is  not  limited  by  the  ability  of 
the  system  to  propagate  a  wave  over  a  large  range,  but  rather  by  the  propagation 
loss  which  increases  with  frequency  [6].  A  target  which  exists  at  some  range  from 
the  radar  site  may  be  illuminated  over  a  frequency  range  of  5  MHz  to  30  MHz; 
however  the  propagation  loss  at  the  upper  frequencies  may  result  in  the  backscatter 
energy  from  the  target  being  on  the  order  of  the  higher-order  sea  clutter  return 
or  galactic  noise  level  therefore  making  detection  and  any  estimation  of  target 
features  impossible.  In  fact,  the  lower  signal  to  noise  ratio  at  higher  frequencies 
may  result  in  poorer  quality  estimates  of  features  even  though  detection  is  still 
possible. 

Interference  with  other  users  in  the  IIF  band  may  also  preclude  the  use  of  some 
frequencies.  During  testing  of  an  experimental  system  in  Florida,  interference  in 
the  8-13  MHz  band  was  a  significant  factor  in  system  operation  [7].  In  addition, 
backscatter  resulting  from  undesired  E  -  region  sky  wave  propagation  periodically 
effected  sytem  performance  at  lower  frequencies  [7].  Therefore  there  are  definite 
limits  to  the  frequency  band  over  which  estimates  of  features  may  be  obtained. 

Thus  once  system  specifications  have  been  set  (type  of  targets,  minimum 
range,  etc.)  the  frequency  range  for  a  surface  wave  radar  may  be  determined. 
Bands  of  6  Mhz  to  8  Mhz  may  expected  to  be  used  [7].  From  here  the  RTI 


system  designer  must  determine  which  features  to  estimate,  at  which  frequencies 
to  make  the  estimates,  and  how  many  total  estimates  are  needed  in  order  to  achieve 
maximum  classification  performance. 

The  frequency  band  over  which  an  HF  sky  wave  radar  can  illuminate  a  target  js 
dependent  on  the  current  ionospheric  conditions.  Therefore  real-time  assessment 
of  the  ionosphere  is  important  in  order  to  determine  the  band.  Assessment  of 
the  frequency  band  may  be  made  through  the  use  of  backscatter  sounding.  A 
backscatter  sounding  consists  of  radiating  a  frequency-modulated  continous-wave 
signal  which  is  swept  in  frequency,  then  receiving  the  backscatter  from  the  earth 
after  it  has  made  two  passes  through  the  ionosphere  [25].  An  oblique  backscatter 
sounding  is  then  produced  by  plotting  the  return  power  level  as  intensity  versus 
time  and  frequency.  An  example  is  given  in  Figure  8. 

A  target  appears  as  a  near  horizontal  line  on  the  resulting  plot,  as  shown  in 
Figure  8.  A  direct  measurement  of  the  frequency  band  over  which  the  target  is 
illuminated  may  then  be  made.  Typical  values  range  on  the  order  of  1  MHz  to 
8  MHz. 

It  must  be  noted  that  real-time  monitoring  of  the  ionospheric  path  in  order  to 
determine  the  largest  frequency  band(s)  over  which  a  target  may  be  observed  are 
not  commonplace.  Other  uses  of  sky  wave  radars  such  as  ocean  waveheight  and 
wind  direction  maps  do  incorporate  real-time  ionospheric  measurements  to  aid  in 
selection  of  frequency  bands  which  are  most  immune  to  multipath  distortion  [21], 
[23],  and  [25]. 

Considerable  research  has  been  and  is  being  conducted  in  the  area  of  real¬ 
time  ionospheric  monitoring.  It  is  apparent  that  these  monitoring  systems  can  be 
of  great  aid  to  an  RTI  system  utilizing  the  IIF  sky  wave  radar  as  its  measurement 
device. 


1 

<1 


i 


a 


i 


3 


F  LAYER 


10  f,  15  fN  20  25 

FREQUENCY  (MHz) 


Figure  8:  Oblique  backscatter  sounding. 


3.3  ABSOLUTE  RCS  ESTIMATION 


The  absolute  radar  cross-section  of  a  target  at  frequency  /  may  be  estimated 


by: 


*(/)  =  PTU ) 


PA\r('TX  ('RX  ^2Fp*CIT  ]  * 

{4ir)3R*Ls 


(3.1) 


where  R  is  the  one  way  path  length  the  wave  travels  to  the  target,  Pay  is  the 
average  transmitted  power,  Gpx  's  ^le  transmitting  antenna  gain,  G  py  ’s  the 
receiving  antenna  gain,  A  is  the  wavelength,  Fp  is  the  one  way  propagation  loss 
factor,  tciT  *s  the  coherent  processing  time,  Pp(f)  is  the  received  power  from 
target.,  and  Lg  is  a  factor  accounting  for  system  losses.  Therefore,  to  directly 
estimate  the  absolute  RCS  of  a  target  at  N  frequencies,  the  above  parameters 
must,  be  estimated  at  each  frequency  used.  System  parameters  such  as  P a \ r ,  G px  i 
G rx i  ^ClTi  f  ^5  are  generally  known  or  may  be  accurately  estimated  at 
each  frequency.  The  power  returned  from  the  target  Pp{J)  at  each  frequency  may 
be  estimated  from  the  appropriate  Doppler  spectra  as  discussed  earlier.  However, 
the  range  ( R )  and  propagation  loss  (Fp)  may  be  much  more  difficult  to  estimate 
depending  on  the  propagation  mode  involved. 

For  a  surface  wave  radar,  the  range  to  the  target  may  easily  be  estimated 
through  time  gating  or  other  signal  processing  strategies.  The  propagation  loss, 
however,  is  a  function  of  the  sea  state  over  the  path  as  well  as  frequency  and 
thus  becomes  more  difficult  to  accurately  estimate;  see  Norton  [6].  A  sky  wave 
radar  poses  a  problem  when  attempting  to  estimate  the  path  length  R.  The  first 
difficulty  results  from  path  length  being  a  function  of  the  ionospheric  structure 
which  is  constantly  changing  with  time.  In  addition,  waves  of  different  frequencies 


may  travel  differenct  paths  to  illuminate  the  same  target  thus  making  H  not  only 
a  function  of  time  but  also  of  frequency  [14]. 

To  circumvent  the  difficulties  of  est  imating  Fp  in  a  surface  wave  radar  and  R 
in  a  skywave  radar,  a  reference  may  be  employed.  Possible  references  include  land 
masses,  ocean  markers,  or  Bragg  lines.  Knowledge  of  the  radar  cross  section  of  the 
reference  then  allows  estimation  of  the  target  RC'S  at  each  frequency.  A  method 
of  calibrating  Bragg  lines  for  use  as  an  absolute  reference  is  given  in  Trisna  [15]. 

3.4  RELATIVE  RCS  ESTIMATION 

Bragg  lines  may  also  be  employed  as  a  relative  reference.  Figure  9  shows  two 
spectra  obtained  at  two  frequencies  f\  and  /2.  Knowing  the  expected  change  in 
the  mean  of  the  Bragg  line  amplitude  with  frequency  (which  is  assumed  to  be  zero 
in  this  example)  the  relative  RCS  estimates  of  the  target  can  be  constructed  as 
shown  in  Eq.  (3.2)  and  Eq.  (3.3). 

<rf(/j)  =  (2°  -  3°)dfl  =  -lOdB  (3.2) 

<4(/2)  =  (15  -  30 )dB  =  -20 dB  (3.3) 

Therefore,  the  RCS  of  the  target  at  f\  is  1  /10th  the  RCS  of  the  reference,  and  the 
RCS  of  the  target  at  fo  is  1/1 00th  the  RCS  of  the  reference.  Thus  the  relative 
change  of  the  target  RCS  has  been  estimated  without  absolute  calibration  of  the 
Bragg  line.  Therefore,  the  use  of  the  relative  RCS  of  the  target  as  a  possible  feature 
for  classification  purposes  should  be  evaluated. 

3.5  RELATIVE  PHASE  ESTIMATION 

Previous  RTI  studies  have  shown  that  intrinsic  target  phase  often  proves  to 
be  a  powerful  feature  when  used  for  classification  [2],  [8].  However,  difficulties  in 
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accurately  estimating  the  path  length  to  the  target  as  in  the  case  of  the  sky  wave 
propagation  may  preclude  the  accurate  estimation  of  intrinsic  target  phase.  There¬ 
fore  an  alternate  method  of  utilizing  target  phase  information  has  b^en  studied. 

Consider  the  phases  of  two  radar  signals  at  frequency  /j  and  /2,  which  are 
returned  from  the  same  target  as  given  by  Eq.  (3.4)  and  Eq.  (3.5). 


2*(Rh)  . 

0/2  =  ■  +  * 


a2 


h 


(3.4) 

(3.5) 


where  R is  the  round  trip  path  length  of  wave  at  frequency  / j,  R ^  is  the  round 
trip  path  length  of  wave  at  frequency  /2,  =  j-,  A2  =  and  9j ]  and  0^  are 

the  intrinsic  target  phases  at  /]  and  /2  respectively. 

By  multiply5 ng  Eq.  (3.4)  and  Eq.  (3.5)  by  Aj  and  A2  respectively  and  then 
subtracting  Eq.  (3.4)  from  Eq.  (3.5),  we  are  left  with: 


A/2  0/2  ~A/i0/i  =  2^Rh  ~  Rf\ )  +  Xh0h  ~  Xh9fl  (3  6) 

where  (Rf2  -  Rfi)  is  the  difference  between  round  trip  path  lengths  at  /2  and 
/ 1  respectively.  Through  careful  selecton  of  A/  —  /2  —  /l  and  using  ionograms 
to  estimate  the  present  structure  of  the  ionosphere,  it  may  be  possible  to  accu¬ 
rately  estimate  {Rf2  -  Rfi  )•  By  multiplying  this  estimate  by  2tt  radians,  and  then 
subtracting  this  result  from  Eq.  (3.6),  we  are  left  with  Eq.  (3.7). 


11  =  Xf2°f2  ~  Xfif,h 


(3.7) 


Further  division  by  A j2  gives  Eq.  (3.8),  which  for  /2  =:  f\  (A/  small)  may  be 
approximated  by  Eq.  (3.9). 


w  =  e 


f  2 


(3.8) 
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w = eh  -  eh 


(3.9) 


Therefore  W  is  a  feature  dependent  on  the  change  of  the  characteristic  phase  of 
the  target  associated  with  /2  and  /j. 

Figure  10  demonstrates  a  physical  interpretation  of  the  W  function  for  a 
simple  scattering  case.  Consider  two  targets,  each  with  single  distinct  scattering 
points  at  f\  and  /2  as  shown.  Each  target  has  an  intrinsic  phase  at  each  frequency 
related  to  the  distance  from  some  reference  to  the  scattering  point  on  the  target 
as  given  in  Eq.  (3.10),  Eq.  (3.11),  Eq.  (3.12)  and  Eq.  (3.13). 
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(3.11) 


(3.12) 


(3.13) 


The  resulting  W  function  for  each  target  is  given  by  Eq.  (3.14)  and  Eq.  (3.15). 
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(3.14) 

(3. 15) 


Note  the  W  function  is  a  measure  of  the  distance  between  the  scattering  for  the 
two  frequencies,  and  not  the  absolute  position  of  the  scattering  point  at  any  one 
frequency. 
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A  multi-frequency  radar  system  illuminating  a  target  with  M  total  adjacent 
frequencies  could  generate  (M  —  1)  W  function  estimates.  These  estimates  could 
then  be  employed  as  additional  features  to  help  identify  the  target.  In  general,  the 
larger  the  A/,  the  more  the  paths  taken  by  the  two  waves  differ,  and  the  harder 
it  is  to  accurately  estimate  (Rf2  —  Rf j  )•  However,  large  A /  values  would  tend  to 
excite  scattering  points  on  the  aircraft  which  are  more  separated  in  distance,  pos¬ 
sibly  leading  to  W  functions  which  contain  more  information  about  the  individual 
targets.  Therefore,  our  ability  to  accurately  estimate  the  change  in  path  lengths 
with  frequency  for  a  given  ionospheric  condition  may  determine  our  maximum  A/. 


CHAPTER  IV 


RADAR  SYSTEM  SIMULATION  AND  EVALUATION 

STUDY 

4.1  INTRODUCTION 

In  order  to  study  the  performance  of  possible  radar  target  classification  sys¬ 
tems  which  utilize  HF  radars  as  a  measurement  device,  a  simulated  system  has 
been  developed.  Over  the  past  several  years  the  ElectroScience  Laboratory  has  de¬ 
veloped  a  large  multi-frequency  data  base  consisting  of  aircraft,  ship,  and  ground 
vehicle  radar  signitures.  This  data  base  has  been  used  to  investigate  several  RTI 
methods  for  the  various  classes  of  targets.  A  data  base  consisting  of  calibrated 
monostatic  radar  returns  of  5  aircraft  was  chosen  for  the  experimental  portion  of 
this  study.  Although  classification  of  the  ships  and  ground  vehicles  is  of  consid¬ 
erable  interest,  the  large  amount  of  previous  work  with  the  aircraft  data  base  [2] 
provides  a  strong  foundation  of  previous  work  for  comparing  results. 

The  ElectroScience  Laboratory  compact  radar  range  facility  was  used  to  create 
the  data  base.  A  basic  block  diagram  of  the  measurement  system  and  a  picture 
of  the  actual  facility  is  given  in  Figures  11  and  12  respectively.  A  more  detailed 
explanation  of  system  may  be  found  in  [3].  In  general,  the  measurement  facility 
simulates  the  illumination  of  a  target  by  a  radar  wave  which  is  approximately  plane 
in  the  region  of  the  target  (locally  plane).  This  would  be  true  of  any  IIF  radar 
system.  The  complex  backscattered  fields  from  the  target  and  other  undesired 
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scatters  are  then  measured  over  a  frequency  band  of  1-12  GHz.  To  obtain  the  high- 
accuracy  range-independent  complex  signiture  of  the  target  ,  a  calibrat  ion  sequence 
described  in  detail  in  Kimball  [17]  was  employed.  Further  detailes  concerning  dat  a 
scaling,  windowing,  and  other  signal  processing  details  are  given  in  Harris  [26], 
and  Kimball  [17].  Five  metallic  scale  models  of  the  Concorde,  DC10,  707,  727,  and 
747  were  used  as  targets.  Silhouettes  and  dimensions  of  the  real  aircraft  are  given 
in  Figures  13  through  17.  Specifications  of  the  5  aircraft  data  base  are  given  in 
Table  1. 

4.2  RADAR  SYSTEM  SIMULATION 

A  flowchart  of  the  Radar  System  Simulation  and  Evaluation  (RSSE)  program 
is  shown  in  Figure  18.  The  program  begins  with  selection  of  N  complex  radar 
signatures  of  the  5  targets  at  frequencies  /],/2v  •  •>/#>  where  the  frequencies  are 
chosen  by  the  operator  to  represent  those  which  may  be  encountered  in  a  typical 
HF  radar  system.  In  addition,  one  azimuth  angle  and  one  polarization  is  selected, 
thus  assuming  these  parameters  have  been  accurately  estimated.  The  target  data 
base  is  then  accessed,  the  data  are  appropriately  scaled  [8],  and  two  sets  of  data 
denoted  the  test  and  catalog  set  are  created.  Each  element  of  the  two  sets  may  be 
represented  by: 

=  Rk{fi)  +  jXk(fi)  (4.1) 


where 


Figure  11:  Compact  range  block  diagram,  (from  Kamis  [2]). 
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Figure  16:  727  silhouette  (from  Kamis  (2]). 


Table  1:  Data  base  specifications  (from  Kamis  [2]). 


Low-Frequency  Data  Base 

Frequency  formatted  data  strings  from  1  to  12  GHz 
Availible  polarizations: 

Transmit  Horizontal,  Receive  Horizontal  (HH) 
Transmit  Vertical,  Receive  Vertical  (VV) 

Transmit  Horizontal,  Receive  Vertical  (VH) 

Availible  Azimuth  angles  @  Elevation  =  0°: 

0°  to  180°  by  10°  and  15°  increments 
Common  aircraft  bandwidth: 

7.6  -  60  MHz  “Scaled” 

Maximum  number  of  usable  frequencies:  209 
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Ak(fi)  —  square  root  of  the  radar  cross  section  of  target  at  frequency 


ft, 

A  corrupted  feature  vector  representing  an  unknown  target  was  constructed 

by  first  generating  2N  total  noise  samples  taken  from  an  independent  Gaussian 

zero-mean  random  process  with  noise  power  given  by  : 

2 

PN  =  101og10(y  )dBsm  (4.4) 

where  cr  is  the  standard  deviation  of  the  process  in  meters  [3].  The  noise  samples 
were  then  added  to  the  N  high  accuracy  complex  radar  signitures  of  a  known 
target  in  the  test,  set  to  create  N  corrupted  complex  data  of  the  form: 


Ru(fi)  =  RT(fi)  +  X(fi) 

(4.5) 

Xu(fi)  =  XT(fi)  +  N'(fi) 

(4.6) 

where  Rj{fi)  is  defined  in  Eq.  (4.2)  for  target  T  in  the  test  set,  Xj( /, )  is  defined 
in  Eq.  (4.3)  for  target  T  in  the  test  set,  and  N(fi),  N'(fj)  are  samples  from  the 
independent  zero-mean  gaussian  noise  process. 

The  resulting  N  complex  corrupted  data  which  may  be  represented  by  : 

Ru(fi)  +  jXu(ft)  =  A(/t)c^'(*),  i  =  1, ...,  N  (4.7) 

constitute  the  unknown  target  feature  vector.  This  corrupted  feature  vector  is 
then  tested  against  the  feature  vectors  of  the  targets  in  the  catalog  set  by  means 
of  a  classifier  function  as  shown.  A  Euclidian  distance  metric  algorithm  commonly 
refferred  to  as  a  Nearest  Neighbor  (NN)  algorithm  was  employed  as  the  classi¬ 
fier  function.  The  NN  algorithm  computes  a  distance  {dv  f.)  as  a  measure  of  the 
distance  between  the  unknown  target  u  and  the  catalog  target  k ,  as  shown  in 


'<m  W,  V '<^<7 V"  X"  yj 


du,k  =  |E(^(/*)-^(/«))‘j  (4-8) 

where  Jri!(ft)  is  the  feature  of  target  k  at  frequency  extracted  from  the  catalog 
set,  and  J-u ( / j )  is  the  feature  of  unknown  target  u  at  frequency  /,  extracted  from 
the  corrupted  feature  vector. 

After  testing  against  the  appropriate  features  of  all  5  possible  cat  alog  targets 
(k  =  1  to  5)  the  unknown  target  is  classified  as  the  catalog  target  which  produced 
the  smallest  distance  measure.  Whether  the  classification  attempt  was  succesful 
or  not  is  known  since  the  identity  of  the  unknown  target  is  known  a-priori.  This 
simulation  process  is  repeated  100  times  for  each  test  target  ,  with  independent 
noise  being  added  each  time  to  create  the  simulated  feature  vector.  Results  of  a 
run  for  a  given  noise  power  are  plotted  as  a  confusion  matrix  as  shown  in  Figure  19. 

The  average  misclassification  percentage  has  a  margin  of  error  since  it  is  the 
result,  of  a  finite,  and  not  an  infinite,  number  of  independent  experiments.  Details 
of  the  error  margin  may  be  found  in  Degroot  [10].  In  general,  the  error  margin  is 
on  the  order  of  2  to  3  percent  for  5  targets  and  100  experiments  [2]. 

4.3  CHANNEL  SELECTION  AND  SIMULATION 


In  genera),  the  ability  to  accurately  estimate  target  features  such  as  abso¬ 
lute  RCS  and  intrinsic  phase  may  be  described  by  defining  channels  with  various 
characteristics  which  affect  the  complex  backscattered  coefficients  of  the  targets. 
This  section  describes  four  channels  with  different  characteristics  which  allow  dif¬ 
ferent  target  features  to  be  estimated.  By  appropriately  processing  the  complex 
backscattered  components  from  the  test  set  for  each  channel  model,  and  alter¬ 
ing  the  complex  backscatter  components  of  the  catalog  set  if  appropriate,  the 
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Figure  19:  Confusion  matrix  (from  Kamis  [2]). 
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resulting  classification  performance  given  the  constraints  of  each  channel  may  be 
synthesized. 

4.3.1  Coherent  Gaussian  Channel  Model  Simulation 

A  channel  which  adds  coherent  independent  zero  mean  Gaussian  noise  to  the 
complex  backscatter  coefficients  of  target  T  is  shown  in  Figure  20.  Note  that 
both  the  amplitude  and  phase  of  the  high-accuracy  backscatter  terms  have  been 
corrupted.  The  resulting  estimates  Ru(fi )  and  X u(fi)  represent  the  corrupted  real 
imaginary  parts  of  the  complex  backscatter  coefficients  of  an  unknown  target  to 
be  classified. 

The  implementation  of  the  Channel  1  model  into  the  RSSE  program  has 
previously  been  described  in  Section  4.2.  An  algorithm  denoted  the  Coherent 
Nearest  Neighbor  Algorithm  given  by: 

'  N 

4,*  =  £(**(/«■)- 

,t=l 

where  A^(/, )  are  the  complex  backscatter  coefficient  of  catalog  target  k, 

and  Ru(fi )»  Xu(fi)  are  the  simulated  real  and  imaginary  part  s  of  complex  backscat¬ 
ter  coefficient  of  unknown  target  u,  is  then  used  to  classify  the  unknown  target  U 
as  being  the  catalog  target  K*  such  that  fc*  is  minimized.  It  is  important  to 
note  that  intrinsic  target  phase,  although  not  directly  estimated  at  the  output  of 
Channel  1,  is  utilized  in  the  estimation  process. 

4.3.2  Non-Coherent  Gaussian  Channel  Model  Simulation 

A  channel  which  corrupts  the  amplitudes  of  the  complex  backscatter  coeffi¬ 
cients  through  coherent  addition  of  independent  zero-mean  Gaussian  noise,  while 
removing  all  intrinsic  phase  information  is  shown  in  Figure  21 .  This  channel  models 


a  radar  system  where  accurate  estimation  of  intrinsic  target  phase  is  not  feasible. 
As  a  result,  only  estimates  of  backscatter  amplitudes  are  available  for  classification. 

The  implementation  of  the  Channel  2  model  into  the  RSSE  program  has 
previously  been  described  in  Section  4.2.  An  algorithm  denoted  the  Non-Coherent 
Nearest  Neighbor  Algorithm  given  by: 

'  JV  15 

du,k=  EiMfi)  ~  Au(ft))2  (4.10) 

L»=i 

where  A^ifi)  is  the  high  accuracy  amplitude  of  complex  backscatter  coefficient  of 
catalog  target  k ,  and  Rv(fi  \  is  the  simulated  amplitude  of  complex  backscatter 
coefficient  of  unknown  target  tt,  is  then  used  to  classify  the  unknown  target  it  as 
being  the  catalog  target  k*  such  that  du  (.*  is  minimized. 

4.3.3  Multiplicitive  Component  Channel  Model  Simulation 

A  channel  which  multiplies  the  amplitudes  of  an  unknown  target’s  backscatter 
coefficients  by  an  unknown  constant  M  in  addition  to  adding  noise  and  removing 
phase  information  is  shown  in  Figure  22.  This  channel  models  a  radar  system 
where  accurate  estimation  of  absolute  RCS  is  not  practical  due  to  the  lack  of  a 
calibrated  reference,  but  estimation  of  relative  RCS  as  described  in  Section  3.4  is 
feasible. 

One  possible  method  of  using  A^(fi)  data  out  of  Channel  3  involves  the 
removal  of  the  unknown  factor  A/.  To  achieve  this,  the  average  power  of  the 
corrupt  vector  out  of  Channel  3  is  calculated  by: 

P*  =  -f  t  Wifi))*  (4.11) 

^  i=  1 

1  N 

pM  =  m2m  D*r(/i)  +  ^'(/())2  +  (AW,)  +  W,))2  (4.12) 

iV  i=i 
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PM  -  m2p 


(4.13) 


where  P  is  the  average  power  of  the  noise  vector  if  M  =  1.  A  new  corrupt  vector 


with  values  given  by: 

AR.n-ilU i> 

ul-M~  vpm 


\A*r(/i)  +  W,))2  +  (*r(/«)  +  iV'(/i))2 


(4.14) 


(4.15) 


1  <  j  <  N. 


(4.16) 


may  then  be  constructed.  This  resulting  vector  now  has  values  independent  of  M 
while  maintaining  the  ratio’s  of  amplitude  estimates  at  different  frequencies,  i.e.: 

W>  -  4./, i  forls,sw'  I£is"-  ,4I6) 

The  Ay(fi)  features  may  now  be  used  for  target  classification  with  a  properly 
chosen  catalog  set. 

The  simulation  of  ut  ilizing  features  passed  through  a  channel  with  an  unknown 
multiplicitive  component  first  required  normalization  of  the  catalog  set.  The  power 
(or  average  RCS)  of  catalog  target  k  is  given  by  Eq.  (4.17). 


h  =  NY.^k{h\ 

e=  1 


(4.17) 


Dividing  each  data  of  catalog  target  k  by  the  square  root  of  P/.  gives  a  new  set  of 


data,  where  the  magnitude  of  each  point  in  the  vector  is  given  by  Eq.  (4.18). 


i  =  \,N 


(4.18) 


The  new  data  A^{ft)  represent  the  relative  radar  cross  section  features  of  catalog 
target  k  which  will  be  used  in  the  classification  routine.  Note  that  the  ratio  of  radar 
cross  section  of  any  two  frequencies  f,  and  fj  for  catalog  target  k  is  preserved  as 
shown  in  Eq.  (4.19). 


I 


Figure  22:  Multiplicative  AWGN  channel  model. 


47 


1  <  t  <  N  ,  1  <  j  <  N 


(1.10) 


(Ak(fi))2  =  Mf(/«))2 
(Mfj))2 

The  normalization  described  in  Eq.  (4.17)  and  Eq.  (4.18)  is  carried  out  on  all 
catalog  targets  (A:  =  1  to  5).  As  a  result,  each  catalog  target  has  an  average  power 


(or  average  RCS)  given  by: 

II 

n'M* 

Q 

to 

(4.20) 

pR_  1  f  (AkUi)V 

Pk  -J»S  \vpJ 

(4.21) 

(4.22) 

P*  =  p/‘  =  I 

(4.23) 

Therefore,  the  overall  average  size  of  each  target  in  the  catalog  has  been  removed. 

The  normalization  is  graphically  depicted  in  Figure  23.  The  complex  backscat- 
tered  amplitudes  at  frequencies  f\  and  for  two  catalog  targets  are  depicted  by 
vectors  V’i  and  Vo.  Vector  V*2  represents  a  relatively  large  target  with  an  average 
RCS  on  the  order  of  1000  sm,  and  V’i  represents  a  smaller  target  with  an  average 
RC'S  on  the  order  of  400  sm.  Note  the  distance  between  the  two  points  in  space 
representing  the  targets  is  on  the  order  of  10  to  15  meters.  After  normalization, 
the  two  targets  are  represented  by  vectors  V7j  and  V^-  Note  the  distance  between 
the  points  in  space  representing  the  two  targets  has  significantly  decreased. 

The  addition  of  coherent  i.i.d.  noise  to  the  complex  backscattered  components 
of  a  test  target  to  form  a  corrupt  test  vector  was  carried  out  as  described  earlier. 
After  the  noise  was  added  the  power  of  the  test  vector  was  calculated  by: 

1  N  0 

P  =  £ (*r(/«)  +  W.:))2  +  (At(/,)  +  nU)Y 


(4.24) 


A' 


The  amplitudes  of  the  simulated  backscatter  components  were  then  normalized 


R...  [(%(/,)  +  N(f ,;))2  +  {XT(fj)  +  7V'(/t))2] j 


AH(fi) 


i  =  1, ...,  Ar(4.25) 


producing  the  normalized  amplitudes  of  the  backscatter  estimates  of  the  unknown 
target.  This  process  would  remove  any  unknown  multiplicitive  component  from 
the  test  vector  as  shown  earlier.  It  is  critical  that  the  noise  be  added  before  the 
normalization  process  is  completed  in  order  to  directly  compare  the  relative  RCS 
to  the  absolute  RCS  as  possible  classification  features.  If  normalization  is  done 
before  noise  is  added,  the  signal  to  noise  ratio  will  be  much  lower  resulting  in 
significantly  poorer  classification  performance.  The  unknown  target  may  then  be 
tested  against  each  of  the  catalog  targets  by  use  of  a  Nearest  Neighbor  Algorithm 
of  the  form: 


(4-26) 


where  A^(ft)  and  are  given  in  Eq.  (4.18)  and  Eq.  (4.25)  respectively,  and 

the  identity  of  the  unknown  target  is  chosen  as  the  catalog  target  k*  for  which  du  ^ 
is  minimized. 

4.3.4  Additive  Component  Channel  Model  Simulation 

Since  the  case  of  using  a  channel  model  with  an  unknown  multiplicitive  factor 
M  is  of  interest,  it  is  natural  to  extend  the  problem  to  the  case  of  an  unknown 
additive  component.  C  as  shown  in  Figure  24.  The  addition  of  the  unknown  com¬ 
ponent  C'2  under  the  radical  sign  as  shown  represents  the  estimation  of  RCS  of  a 
target  at  2  or  mere  frequencies,  where  only  the  difference  in  RCS  is  known.  For 
example,  the  difference  in  RCS  of  a  target  at  /]  and  fo  may  be  known  to  be  500  sm 


Vi 

& 


>/;• 

S: 

s* 


w 


but  whether  the  RC’S  at  f\  and  fi  is  10  sm  and  510  sni  respectively,  or  1000  sm 
and  1500  sm  is  not  known.  Although  this  problem  may  or  may  not  arise,  it  is  of 


>*?(/.)  =  [(flr(/,)  +  iV(/,))2  +  (A-r(/,)  +  A'(/.))2 

+  C2  -  [P„orm  +  (C2  -  P)}]'  (4.30) 

A?(/,)  =  [(ftr(/,)  +  N(fi))2  +  (XT (f,)+ N'(f,))2  +  [P  -  Pnom,]]’  (4.31) 

giving  a  new  vector  independent  of  C,  while  maintaining  the  difference  in  power 
(RCS)  between  data  at  any  two  frequencies. 

Features  constructed  from  a  channel  with  an  unknown  constant  additive  com¬ 
ponent  are  referred  to  as  relative  difference  RCS  (RDRCS)  features.  The  features 
were  implemented  into  the  RSSE  program  by  reconstructing  each  target  in  the 
catalog  set  so  its  new  average  power  '  is  the  same  for  all  targets  while  the  dif¬ 
ference  between  RCS  at  a  given  frequency  and  the  average  RCS  is  maintained.  The 
same  reconstruction  of  the  corrupted  test  vector  removes  any  unknow'n  additive 
component  from  the  estimates.  As  a  result,  each  target  in  the  catalog  set  is  given 
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Figure  24:  Additive  AWGN  channel  model. 


the  same  average  RCS,  with  the  difference  between  the  RCS  at  each  frequency  and 
the  average  RCS  being  maintained. 

The  normalization  of  the  catalog  set  in  the  RSSE  program  was  accomplished 
as  described  above.  Since  the  amplitude  of  the  complex  1;:  Waiter  coefficients 
in  the  catalog  set  is  actually  the  square  root  of  the  radar  cross  section  at  each 
frequency,  the  amplitudes  were  reconstructed  using: 

A^fi)  =  [(Mfi)f  +  (Pnorm  ~  Pk)]'*  .  (4.32) 

where 

1  N 

h  =  ^E(Ak(fi))2-  (4.33) 

JV  1=1 

An  example  of  the  normalization  effect  on  two  catalog  targets  is  given  in 
Figure  25.  Consider  target  1  in  the  cat  alog  set  with  RCS  values  500sm,  400sm,  and 
300sm,  and  target  2  with  RCS  values  of  600sm,  lOOOsm,  and  1400sm  at  frequencies 
/ 1,  /2  and  as  shown.  The  average  RCS  over  the  three  frequencies  is  400sm  and 
lOOOsm  for  targets  1  and  2  respectively.  After  the  normalization,  the  average  RCS 
of  both  targets  is  IGOOsm  ( Pnorm  =  1600sm)  as  shown.  As  a  result,  absolute 
target  size  is  lost;  however,  the  difference  between  the  RCS  at  any  frequency  and 
the  average  RCS  is  maintained. 

Coherent  noise  was  added  as  always  to  the  complex  backscatter  coefficients  of 
target  T  of  the  test  set  to  create  the  noisy  feature  vector.  The  average  power  of 
the  noisy  feature  vector  was  then  calculated  as: 

1  N 

PV  =  v  E(4(/i))2-  (4-34) 

1  i=l 

The  amplitudes  of  the  noisy  feature  vector  were  then  reconstructed  as: 

AS(fi)  —  \i.Av{fi))“  +  (Pnorm  ~  PU )]  "  (4.35) 


FREQUENCY 


X  TARGET  I  ABSOLUTE  RCS 
H  TARGET  I  NORMALIZED  RCS 

O  TARGET  2  ABSOLUTE  RCS 
B  TARGET  2  NORMALIZED  RCS 


Figure  25:  Normalization  of  two  target 


These  data  represent  the  estimates  of  the  amplitudes  of  the  complex  coefficients 
passed  through  channel  4  after  the  additive  unknown  component  has  been  removed. 


As  a  result,  the  Nearest  Neighbor  Algorithm  given  by: 

l 


(4.36) 


where  Af?(ft)  and  A^(ft)  are  given  by  Eq.  (4.32)  and  Eq.  (4.35)  respectively,  was 
used  to  classify  the  unknown  target  u  as  the  catalog  target  k*  for  the  k *  such  that 
dujf*  is  minimized. 

It  is  important  to  note  that  although  the  power  of  the  catalog  vectors  and 
corrupted  vector  is  altered,  the  signal-to-noise  ratio  of  the  init  ial  corrupt  ed  vector  is 
not  changed,  and  thus  classification  performance  will  not  be  function  of  the  power 
level  chosen  for  normalization.  Simulations  were  repeated  using  normalization 
power  levels  of  4000  sm,  6000  sm,  and  8000  sm  to  insure  this  result  . 


4.4  ASSESSMENT  OF  MISCLASSIFICATION  CURVES 


A  typical  set  of  misclassification  curves  is  shown  in  Figure  26.  Performance  of 
a  given  curve  with  respect  to  another  may  be  related  by  the  relative  noise  power 
immunity.  The  noise  power  immunity  of  a  curve  is  determined  by  the  difference 
between  the  noise  power  of  the  curve  and  the  reference  curve  evaluated  at  a  con¬ 
stant  misclassification  level  (generally  between  5  and  40  percent).  The  noise  power 
immunity  is  an  approximate  measure  of  the  additional  noise  power  which  must  cor¬ 
rupt  the  according  data  to  produce  a  curve  equivalent  to  the  reference.  A  positive 
noise  power  immunity  indicates  better  classification  performance  with  respect  to 
the  reference,  while  negative  noise  power  immunity  indicates  the  opposite. 
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CHAPTER  V 


RSSE  RESULTS  FOR  AWGN  CHANNELS 

5.1  Frequency  Band  Study 

The  frequency  band  over  which  a  target  at  a  given  range  may  be  illuminated 
is  limited.  As  a  result,  only  estimated  features  from  the  limited  band  are  available 
for  classification.  Thus  a  study  has  been  conducted  to  find  the  effect  of  frequency 
band  size  on  classification  performance  for  the  four  channels  described  in  Sections 

4.3.1  through  4.3.4. 

Figure  26  shows  the  change  in  performance  as  the  frequency  band  of  channel  1 
is  decreased.  Note  the  overall  performance  is  not  significantly  effected.  Figure  27 
demons4 rates  the  greater  loss  of  performance  with  decreasing  frequency  band  when 
channel  2  used.  A  significant  loss  of  noise  power  immunity  of  approximate!  8  dB-sm 
is  encountered  when  the  frequency  band  decreased  from  4  MHz  to  2  MHz. 

Similarly,  Figure  28  demonstrates  the  even  greater  decrease  in  classification 
performance  encountered  by  decreasing  the  frequency  band  when  channel  3  intro¬ 
duces  an  unknown  multiplicitive  component.  Losses  on  the  order  of  5  and  12  dB-sm 
noise  power  for  a  given  misclassification  percentage  appear  as  the  frequency  band 
is  decreased  from  8  MHz  to  4  MHz,  and  then  4  MHz  to  2  MIIz.  Figure  29,  which 
was  generated  utilizing  the  unknown  additive  component  channel,  expresses  the 
most  significant  decrease  in  performance  as  the  band  is  decreased  from  8  MHz  to 


Figures  26  through  29  demonstrate  some  qualilies  unique  to  each  of  the  four 
channels.  Channel  1  demonstrated  a  very  small  change  in  performance  compared 
with  the  three  other  channels  as  the  frequency  band  was  significantly,  decreased. 
Both  channel  2  and  channel  3  showed  a  much  larger  degredation  in  performance 
when  the  band  was  reduced  to  2  MHz  as  compared  to  the  reduction  to  4  MHz. 
This  implies  a  frequency  band  of  not  less  than  4  MHz  would  be  highly  desirable 
given  a  channel  which  does  not  allow  the  estimation  of  phase  information.  This 
also  implies  that  a  doubling  of  the  frequency  band  from  4  MHz  to  8  MHz  does  not 
provide  the  same  degree  of  improvement  in  classification  performance.  Therefore 
the  cost  of  having  the  ability  to  obtain  features  over  an  8  MHz  frequency  band,  such 
as  a  second  radar  site,  may  not  be  warranted.  Conversely,  channel  4  demonstrated 
significant  performance  degredation  in  both  cases  of  band  reduction.  Therefore, 
given  a  channel  with  an  unknown  additive  component,  a  significant  increase  in 
performance  may  be  obtained  by  extending  the  frequency  band  beyond  4  MHz. 
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Parameter  Under  Observation:  Frequency  Band  Size 

Channel  Model:  Non- Coherent  AWCN  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  Varies 

Number  of  Frequencies:  10 


_ 8.0  MHz  to  10.0  MHz  Band 

-  -  -  -  8.0  MHz  to  12.0  MHz  Band 

—  -  —  8.0  MHz  to  1G.0  MHz  Band 


Figure  27:  Noil-Coherent  channel  band  size  study. 
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Parameter  Under  Observation:  Frequency  Band  Size 

Channel  Model:  Multiplicitive  Component  AVVGN  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  Varies 

Number  of  Frequencies:  10 


_ 8.0  MHz  to  10.0  MIlz  Band 

- 8.0  MHz  to  12.0  MHz  Band 

- 8.0  MHz  to  16.0  MHz  Band 


Figure  28:  Multiplicitive  component  channel  band  size  study- 
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PERCENT  MISCLASSFICATION 
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Parameter  Under  Observation:  Frequency  Band  Size 

Channel  Model:  Additive  Component  AWCN  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MIIz 

Stop  Frequency:  Varies 

Number  of  Frequencies:  10 
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NOISE  POWER  (DB-SM) 


_ 8.0  MHz  to  10.0  MHz  Band 

-  -  -  -  8.0  MHz  to  12.0  MHz  Band 

—  -  —  8.0  MHz  to  16.0  MHz  Band 


Figure  29:  Additive  component  channel  band  size  study. 
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5.2  Number  of  Features  Study 


A  study  was  conducted  to  determine  the  relationship  between  the  number  of 
features  used  and  the  classification  performance  for  the  four  channels.  Previous 
research  [2]  concluded  that  the  minimum  frequency  sampling  interval  which  will 
satisfy  Shannon’s  sampling  criteria  is  given  by: 

Afisk  '■5I) 


where  c  is  the  speed  of  light  and  L  is  the  length  of  the  target.  Therefore  given  an 
average  target  length  of  50  meters,  one  finds  approximately  2,  3,  and  5  frequency 
samples  are  required  to  meet  the  above  criteria  for  band  widths  of  2  MHz,  4  MHz, 
and  8  MHz  respectively.  These  values  were  used  as  minimum  guidelines  to  be  met 
for  each  corresponding  frequency  band.  Since  Kamis  [2]  examined  the  effect  of  the 
number  of  features  at.  large  frequency  bands  (Fian(j  >  8  MHz)  for  channels  1  and 
2,  these  bands  for  these  channels  were  not  included  in  the  study. 

Figure  30  shows  a  significant  increase  in  noise  immunity  of  approximately 
5  dB-sm  as  the  number  of  features  through  channel  1  are  increased  from  3  to  5 
to  10  over  the  2  MHz  band.  The  increase  in  noise  immunity  of  approximentally  5 
dB-sm  demonstrates  that  classification  performance  can  be  significantly  improved. 
Figure  31  exhibits  the  smaller  noise  immunity  of  approximately  3dB-sm  associated 
with  channel  2  as  the  number  of  features  is  increased  from  3  to  10.  Figure  32 
demontrates  an  almost  constant  increase  of  noise  immunity  of  3  dB-sm  to  5  dB-sm 
as  the  number  of  features  is  increased  from  3  to  5  to  10  to  20  in  the  channel  with  an 
unknown  multipicitive  component  (channel  3).  Figure  33  demonstrates  a  similar 
increase  in  performance  for  the  unknown  additive  component  channel. 

Figures  30  through  33  illustrate  that  classification  performance  may  be  im¬ 
proved  by  increasing  the  number  of  features  when  limited  to  a  2  MHz  frequency 


band,  irregardless  of  the  channel.  Channel  2  demonstrated  the  smallest  increase  in 
performance  of  the  three,  however,  the  increase  is  much  greater  than  that  observed 
by  Kamis  [2]  when  increasing  the  number  of  features  over  larger  frequency  bands 
of  8  MHz  and  22  MHz. 

Figures  34  and  35  show'  the  relationship  between  the  number  of  features  and 
classification  performance  over  4  MHz  and  8  MHz  frequency  bands  for  channel  3. 
Both  figures  indicate  an  additional  noise  immunity  of  approximated  3  to  5  dB  as 
the  number  of  features  is  increased  from  5  to  10.  A  smaller  increase  in  perfor¬ 
mance  is  obtained  by  increasing  the  number  of  features  from  10  to  20  and  10  to 
15  as  shown  in  Figures  34  and  35  respectively.  Therefore,  as  frequency  band  is 
increased,  channel  3  demontrates  a  decrease  in  the  amount  of  improvement  in  clas¬ 
sification  performance  obtained  by  using  an  increased  amount  of  features.  Note 
that  channel  2  displays  this  same  property.Figure  36  demonstrates  that  the  addi¬ 
tive  component  channel  exhibits  the  same  properties  as  channel  2  and  channel  3 
concerning  frequency  band,  number  of  samples,  and  classification  performance. 

Therfore,  Figures  30  through  36  indicate  that  if  one  is  limited  to  a  small 
frequency  band,  it  may  be  worth  the  cost  of  the  additional  processing  time  to 
obtain  estimates  at  10  or  20  frequencies  regardless  of  the  channel.  However,  if  a 
larger  frequency  band  is  available,  the  additional  processing  cost  may  outweigh 
the  slight  increase  in  classification  performance  obtained  through  the  additional 


Parameter  Under  Observation:  Number  of  Features 

Channel  Model:  Coherent  AWGN  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MIIz 

Number  of  Frequencies:  Varies 


Parameter  Under  Observation:  Number  of  features 

Channel  Model:  Non-Cohcrent  AYVGN  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  Varies 


_  3  Features 

-  -  -  -  5  Features 

—  -  —  10  Features 


Figure  31:  Non-Coherent  channel  number  of  features  study  (8.0-10.0 
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Parameter  Under  Observation:  Number  of  Features 

Channel  Model:  Multiplicitive  Component  AWGN  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  Varies 
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Figure  32:  Multiplicitive  component  channel  number  of  features  study 

(8.0-10.0  MHz). 
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Parameter  Under  Observation:  Number  of  Features 

Channel  Model:  Multiplicitive  Component  AWGN  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  Varies 
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Figure  34:  Multiplicitive  component  rhannel  number  of  features  study 

(8.0-12.0  MHz). 
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5.3  Comparison  of  AWGN  Channels  Performance 


A  direct  comparison  oi  the  lour  channels  lias  been  made  in  terms  of  clas¬ 
sification  >erformance.  Figure  37  demonstrates  the  dependence  of  classification 
performance  on  the  channel  given  a  2  MHz  frequency  band.  A  significant  perfor¬ 
mance  degredation  of  at  least  15  dBsm  is  experienced  for  channel  2  with  respect  to 
channel  1.  The  classification  performance  of  channels  3  and  4  is  further  degraded 
as  a  result  of  target  information  lost  in  the  normalizat  ion  processes  used  to  remove 
the  unknown  (multiplicitive  or  additive)  components. 

Figure  38  shows  the  resulting  performance  of  the  same  channels  over  a  wider 

4  MHz  frequency  band.  The  performance  of  channels  2,  3,  and  4  improves  sig¬ 
nificantly  with  respect  to  the  performance  of  channel  1.  At  lower  noise  powers, 
channel  3  is  found  to  outperform  channel  4  by  a  surprisingly  wide  margin  of  ap¬ 
proximately  2  dBsm. Figure  39  compares  the  performance  of  the  4  channels  over  an 
8  MHz  bandwidth.  The  additive  component  channel  now  outperforms  the  multi¬ 
plicitive  channel.  Note  that  the  performance  of  channels  2.  3,  and  4  all  fall  within 

5  dBsm  of  each  other. 

Figures  37,  38,  and  39  together  indicate  that  channels  with  unknown  mul¬ 
tiplicitive  or  additive  components  which  must  be  removed  will  always  degrade 
performance.  This  is  understandable,  since  information  concerning  target  size  is 
lost  in  the  normalization  process  used  to  remove  the  unknown  (multiplicitive  or  ad¬ 
ditive)  component. More  importantly,  these  figures  indicate  that  the  performance 
of  a  channel  with  an  unknown  additive  component  (channel  4)  degrades  at  a  faster 
rate  with  a  narrowing  frequency  band  than  a  channel  with  an  unknown  multi¬ 
plicitive  component  (channel  3).  Figures  37  and  38  clearly  indicate  that  knowing 
the  difference  in  absolute  RCS  between  frequencies  (channel  4)  will  not  provide 
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any  better  performance  over  knowing  only  the  relative  change  in  ROS  (channel  3). 
Finally,  the  three  figures  indicate  that  without  phase  information  classification 
performance  decreases  dramatically  with  a  decrease  in  frequency  band. 

Figures  40  and  41  show  the  classification  performance  of  the  four  channels 
using  5  and  15  features  respectively  in  the  8  MHz  frequency  band.  Again,  chan¬ 
nel  1,  which  utilizes  phase  information,  performs  significantly  better.  Note  that 
the  relative  performance  of  each  channel  with  respect  to  the  other  three  remains 
almost  constant  in  both  figures.  Therefore  all  three  channels  tend  to  suffer  the 
same  loss  of  performance  with  the  lowering  of  the  number  of  frequencies  in  this 
large  8  MHz  frequency  band. 

Figures  42  displays  the  resulting  poor  performance  of  channels  3  and  4  when 
the  number  of  features  is  limited  to  three  and  the  frequency  band  (2  MHz)  is 
limited.  The  9  dB-sm  increase  in  noise  power  immunity  of  channel  2  with  respect 
to  channel  3  at  the  20%  misclassification  level  is  one  of  the  largest  changes  in 
performance  found  between  the  two  channels. 

Figure  43  demonstrates  the  relatively  large  improvement  in  classification  per¬ 
formance  of  channel  3  with  respect  to  channel  2,  especially  at  noise  power  levels  of 
less  than  5  dB-sm.  This  increase  in  relative  classification  performance  demostrates 
that  the  performance  of  channel  3  degrades  faster  than  the  performance  of  chan¬ 
nel  2  when  the  number  of  features  is  decreased.  This  same  property  was  noted  in 
the  number  of  feat  ures  study.  Finally,  Figure  44  demonst  rates  an  even  smaller  de¬ 
crease  in  classification  performance  for  channel  3,  further  backing  the  performance 
loss  —  number  of  frequency  relationship  for  channel  3  at  small  band  widths. 

Figures  45  and  46  compare  the  channels  over  a  different  band  from  20  MHz  to 
24  MHz  and  20  MHz  to  28  MHz.  The  channels  demonstrate  the  same  characteris¬ 
tics  over  these  bands  as  over  the  8  MHz  to  12  MHz  and  8  MHz  to  16  MHz  bands. 
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Therefore  the  channel  characteristics  which  have  been  noted  do  not  appear  to  be 
dependent  on  the  exact  band  in  the  HF  spectrum  used,  but  only  depend  on  the 
size  of  the  band. 
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Parameter  Untler  Observation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 


_ Coherent  Channel 

—  -  -  -  Non-Coherent  Channel 

—  -  —  Multiplicitive  Component  Channel 

—  —  —  Additive  Component  Channel 


Figure  37:  Channel  study  over  2  MHz  band. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 
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—  -  -  -  Noil-Coherent  Channel 
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Figure  38:  Channel  study  over  4  MHz  band. 
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Parameter  Under  Observation: 
Channel  Model:  Varies 
Azimuth  angle:  45 
Start  Frequency:  8.0  MIIz 
Stop  Frequency:  16.0  MHz 
Number  of  Frequencies:  10 


Channel 
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Coherent  Channel 
Non-Cohereut  Channel 
Multiplicitive  Component  Channel 
Additive  Component  Channel 


Figure  39:  Channel  study  over  8  MHz  band. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  16.0  MHz 

Number  of  Frequencies:  5 


NOISE  POWER  (DB-SM) 

_ Coherent  Channel 

-  -  Noil-Coherent  Channel 
—  Multiplicitive  Component  Channel 


Figure  40:  Channel  study  for  r>  frequency,  8  MHz  band. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Varies 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  3 
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Figure  42:  Channel  study  for  3  frequency,  2  MHz  band. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 
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Figure  44:  Channel  study  for  10  frequency,  2  MHz  band. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  45 

Start  Frequency:  20.0  MHz 

Stop  Frequency:  24.0  MHz 

Number  of  Frequencies:  10 
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Figure  45:  Channel  study  for  upper  4  MHz  band 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  45 

Start  Frequency:  20.0  MIIz 

Stop  Frequency:  28.0  MHz 

Number  of  Frequencies:  10 
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Figure  40:  Channel  study  for  upper  (8  MHz  band. 
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CHAPTER  VI 


RELATIVE  PHASE  CHANNEL  MODEL 

6.1  Feature  Implementation 

The  development  of  the  IV  feature  in  Section  3.5  was  intended  to  provide 
information  related  to  the  intrinsic  phase  of  the  target  at  times  when  the  accurate 
estimation  of  absolute  path  length  is  difficult  to  achieve.  It  has  previously  been 
shown  that  given  certain  channel  conditions  one  may  create  ( Ar  —  1 )  target  features 
from  the  difference  of  the  intrinsic  phases  at  adjacent  frequencies.  These  (N  —  1) 
features,  coupled  with  N  estimates  of  the  RCS  of  the  target  at  each  frequency, 
provide  a  total  (2 N  —  1)  features  of  the  target  which  may  be  used  for  classification 
purposes.  The  (N  —  1)  W  features  contain  additional  information  about  the  target 
provided  they  are  not  too  severely  affected  by  noise.  If  used  correctly  in  a  clas¬ 
sification  algorithm  along  with  the  amplitude  features,  these  IV  features  should 
improve  classification  performance  over  that  of  using  the  amplitude  features  only. 

A  variation  of  the  nearest-neighbor  algorithm  was  developed  to  incorporate 
the  W  features  in  the  classification  process.  Previous  experience  has  shown  the 
nearest-neighbor  algorithms  can  provide  powerful  classification  performance  results 
[2],  and  therefore  have  been  used  as  the  starting  point.  In  addition,  the  use  of  a 
distance  metric  makes  intuitive  sense  in  a  problem  of  this  nature.  Consider  two 
catalog  targets  (k  =  1,  2)  whose  backscatter  amplitudes  ( \/ RC'S)  measured  at  / j, 
/2,  and  fo  are  plotted  in  3  dimensional  space  as  in  Figure  47.  These  two  targets 


have  two  corresponding  intrinsic  II'  values  each  given  by  Eq.  (6.1 )  and  Eq.  (6.2) 
which  may  be  plotted  in  2-dimensional  space  as  shown  in  Figure  48. 


h(h), 

k  =  1,2 

(6.1 

@k(f2  )i 

A-  =  1,2 

(6.2 

where  0j.(/t)  =  intrinsic  phase  of  target  I\  at  frequency 

Next  consider  estimating  the  corrupted  amplitude  and  corrupted  IF  values  of 
an  unknown  target  u  at  the  corresponding  frequencies  / j,  /2,  and  and  plotting 
the  estimates  in  Figure  47  and  Figure  48  as  shown.  Use  of  a  Euclidian  distance 
metric  provides  four  distinct  measures  as  given  in  Eq.  (6.3),  Eq.  (6.4),  Eq.  (6.5), 
and  Eq.  (6.6). 

t 

3  5 

<1=  (fi)-AMi))2  (6.3) 

i=l 


<2=  yZ(A2ifi)-Au[fi)y 


<1=  : Lm(fi)-wu(fi)r 

.1=1 


<2=  DH'2(/i)  -»«(/<))■ 


A  final  distance  between  the  unknown  target  and  the  two  cat  alog  targets  may 
then  be  constructed  as: 


jA  +  W  _  jA  ,  rtjlf 

«u,l  -  rfu,l  +  PWdu,l 


(meters) 


^u,2  W  ~  ^u,2  +  PQd™. 2 


(meters) 


(6.8) 


where  the  p  term  provides  a  weighting  of  the  W  metric,  while  the  Q  term  provides 
a  unit  conversion  from  an  angle  measure  to  a  length  measure. 

The  Q  term,  labeled  the  catalog  normalization  factor,  is  a  measure  of  the 
separation  of  the  targets  in  the  catalog  set.  In  Figure  47  ,  the  average  separation 
of  the  two  catalog  targets  in  terms  of  amplitude  features  may  be  calculated  by 
Eq.  (6.9). 

l 

3 

Likewise,  the  average  separation  of  the  catalog  targets  in  terms  of  IF  features  may 
be  calculated  by  Eq.  (6.10). 

2  12 

r>C  =  £(H'2(/.)-H',</i))2  (degrees)  (6.10) 

.‘=1 

The  ratio  of  these  two  terms,  as  given  by  Eq.  (6.11),  provides  the  catalog 
normalization  factor. 


Qc  =  pw 
uc 


(6.11) 


The  catalog  normalization  factor  may  be  viewed  in  the  following  manner.  Fig¬ 
ures  47  and  48  show  the  same  two  catalog  targets  in  terms  of  two  different  met¬ 
rics  separated  by  two  different  average  distances  as  calculated  in  Eq.  (6.9)  and 
Eq.  (6.10).  By  multiplying  the  IF  features  of  the  catalog  set  by  Q  (rescaling  the 
phase  feature  space),  one  can  find  the  new  average  separation  of  the  catalog  targets 
in  terms  of  the  IF  features  as: 


d'c  =  Q2E(»2(/,)-vri(/i))2 

»=i 

/>jr  =  q  -  Djy 


(6.13) 


(6.14) 


/>jr  =  d& 


Therefore  multiplying  the  catalog  IT  features  by  Q  lias  the  effect  of  convert¬ 
ing  the  IT  feature  space  to  an  amplitude  feature  space  by  using  average  target 
separation  as  the  conversion  factor.  Likewise,  the  rescaling  of  d by  Q  as  shown 
in  Eq.  (6.7)  and  Eq.  (6.8)  has  the  same  resulting  effect. 

The  value  of  p,  the  IT  feature  weight,  reflects  the  emphasis  placed  on  the  IT 
features  in  the  classification  process.  For  p=0,  the  IT  features  are  ignored,  and  the 
classification  decision  is  solely  determined  by  the  amplitude  distance.  For  p  >  1, 
the  normalized  IT  features  are  more  heavily  weighted  than  the  amplitude  features. 

In  order  to  determine  a  proper  value  of  p,  one  must  consider  the  radar  system 
from  which  the  estimates  of  amplitude  and  phase  are  made.  In  general,  the  de¬ 
gree  of  corruption  of  the  IT  features,  with  respect  to  the  amplitude  features,  will 
determine  the  emphasis  or  de-emphasis  of  the  normalized  IT-distance  measure. 

6.2  Channel  Simulation 


The  relative  phase  channel  model  was  simulated  as  follows.  The  IT  values 
of  the  catalog  set  were  first  created  using  the  complex  backscattered  coefficients 
previously  loaded  in  the  catalog  set.  The  normalization  factor  of  the  catalog  set 
was  then  created  by: 


Qcatalog  ~ 


sLiSiU-H  Tfii 

eL,  £U+,  [s"T  cm/,)  -  iw.))s] 


16.10) 
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where  Aa(fi),  Ai(fi)  are  the  barkscatter  amplitudes  of  catalog  targets  a  and  6, 
and  Hra(/j),  H^/j)  are  the  IT  values  of  catalog  targets  a  and  b.  A  corrupt  test  set 
vector  was  formed  through  the  addition  of  coherent  zero  mean  Gaussian  noise  to 
a  high  accuracy  test  set  vector  as  described  earlier.  The  (AT  -  1)  noisy  IT  features 
were  then  created  from  these  coherent  backscatter  amplitudes  by: 


Ur(/,+.)/  *;+>  V«T(/,)/ 


i  =  1,7V  -  1(6.17) 


where 


AX/,)  =  AX/,)  +  AT',  and 

RT(fi)  =  RTift)  +  N", 


given  Rx(fi  )>  ^T(fi)  are  the  high  accuracy  complex  backscatter  coefficients  of 
test  target  T  at  frequenct  /j  (  or  /,+  i),  and  N1 ,  N"  are  i.i.d.  Gaussian  noise 
samples.  Note  that  the  distribution  if  Hr„(/i)  is  going  to  be  dependent  on  Ax(fi) 
and  the  distribution  of  Au{fi).  This  may  be  visualized  in  Figure  49.  This  graphic 
description  of  the  addition  of  coherent  noise  to  a  large  and  small  target  shows  the 
strong  dependence  of  standard  deviation  of  liru(/j)  on  the  uncorrupted  amplitude 
Ajifi )  of  the  data  point.  Therefore  targets  with  large  RCS  values  are  biased 
toward  having  W  features  that  are  less  corrupted  for  a  given  noise  power. 

An  understanding  of  the  relationships  between  the  distributions  of  11  „(/,)  and 
Au(/i)  for  given  values  of  Ax(fi)  and  levels  of  noise  power  is  important  in  order 
to  determine  sensible  values  of  p.  To  accomplish  this,  a  study  was  conducted  by 
first  selecting  targets  with  a  wide  range  of  backscatter  amplitudes.  The  standard 
deviation  of  the  IF  features  were  estimated  as  described  in  Degroot  [10]  by: 


where  Wj(fi)  is  the  high  accuracy  If  feature  ami  lEj .(/j)  is  the  rth  independent 
corrupt  W  feature.  The  min  statement  specifies  the  minimum  angle  between  the 
values  to  eliminate  branch  cut  effects.  In  a  similar  manner,  the  standard  deviation 
of  the  corrupted  amplitude  features  was  calculated  by  Eq.  (6.19). 


(6.19) 


aA(fi)  =  ™ 


Estimates  °f  )  versus  Mfi)  at  noise  power  levels  of  10  dB-sm  and  20  dB- 
sm  are  shown  in  Figure  50.  Note  that  is  higlily  dependent  on  both  A(f, ) 

and  the  noise  power.  As  a  result,  one  may  desire  to  weight  the  distance  between  an 
unknown  IE  feature  at  frequency  /„  and  the  catalog  IE  feature  at  fn  dependent 
on  the  amplitude  feature  Au  at  that  frequency.  The  standard  deviation  of  the 
amplitude  feature  (<r  t  )  was  found  to  be  almost  independent  of  A(fi)  and  only 


dependent  on  noise  power.  Estimated  values  were  <r 


0.3  sm  and  1.0  sm  at 


noise  power  values  of  10  dBsm  and  20  dBsm,  respectively.  The  previous  study 
led  to  the  development  of  two  IE  feature  distance  measurements.  The  standard 
amplitude  distance  measurement  given  by: 


(6.20) 


was  used  to  measure  the  distance  between  cat  alog  target  k  and  unknown  target  xi 
in  terms  of  amplitude  features.  One  IE  distance  measurement  was  given  by: 

1  E'('n(/i)-H'„(/,))T  (6.21) 

t=l 

where  Hrj(.(/,)  and  Wu{fi)  represent  the  catalog  and  corrupted  IE  features  respec¬ 
tively.  A  second  IE  distance  measurement  was  given  by: 


2 =  £  R(A{fi))  (Wfc(/«)  -  VE„(/,)): 


(6.22) 
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STANDARD  DEVIATION  ESTIMATES  (DEGREES) 


where  R{A{fi))  is  a  function  of  the  corrupted  amplitude  of  t  lie  test  vector.  R(A(fi )) 
is  a  weight  which  approaches  one  for  large  values  of  ,4(/,)  (<r^^  )  sma^)  an(l  aP~ 
proaches  0  for  small  values  of  A{fj)  large).  Note  that  the  amplitudes  of 

the  corrupted  vector  must  be  used  and  not  the  amplitudes  of  the  catalog  presently 
being  tested  agianst.  Use  of  the  uncorrupted  catalog  amplitudes  would  bias  2d]\. 
to  be  small  for  small  targets  and  large  for  large  targets,  therefore  biasing  the 
classification  decision. 


6.3  Experimental  Results 


Classification  runs  were  made  using  Eq.  (6.21)  for  \  to  determine  the 
effects  of  p  on  the  classification  algorithm.  Figures  51,  52,  and  53  show  the  results 
of  nine  tests  run  over  three  separate  frequency  bands.  Figure  51  demonstrates 
an  average  improvement  of  6  dB-sm  of  noise  power  immunity  using  the  weight  of 
p  ~  1.0  over  the  small  2  MHz  band.  However,  Figures  52  and  53  demonstrate  a 
loss  in  performance  as  the  k  distance  weight  is  increased.  These  results  indicate 
that  the  id^j  f,  distance  measure  is  less  than  optimal  since  it  weights  all  IT  features 
equally,  irregardless  of  the  amount  of  corruption  of  the  IT  feat  ure. 

A  second  set  of  tests  were  performed  to  investigate  possible  functions  for 
R(A(fi))  given  in  Eq.  (6.22).  In  this  set  of  tests  the  weighting  was  set  at  p  —  0.5. 
The  two  functons  chosen  for  R(A(fi))  were: 

Mh) 


R(A(fi))  = 


36.0  ’ 


(6.23) 


and 


R(A{fi))  = 


Ml 


36 


»• 


(6.24) 
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since  a  value  of  A(fi)  —  36  meters  corresponds  to  a  RC'S  of  approximately  1300 
square  meters,  or  31  dB-sm,  corresponding  approximately  to  the  largest  RC'S  value 
in  the  5  target  library  at  HF  frequencies. 

Figure  54  demonstrates  improvement  in  classification  performance  for  the  rel¬ 
ative  phase  channel  over  the  noil-coherent  channel,  however  the  increase  in  perfor¬ 
mance  is  less  than  that  shown  in  Figure  51.  Figures  55  and  50  show  little  or  no 
improvement  in  classification  performance  of  the  relative  phase  channel  over  the 
non-coherent  channel  when  utilizing  the  2^v  k  measure.  Therefore  our  choices  of 
R(A (ft)}  are  far  from  ideal. 

Overall  this  study  has  indicated  that  for  small  frequency  bands,  the  relative 
phase  channel  described  may  improve  performance  over  the  non-coherent  channel 
if  the  relative  phase  estimates  are  used  as  described.  However,  at  larger  frequency 
bands,  the  use  of  relative  phase  information  proposed  appears  only  to  deteriorate 
performance.  This  may  be  due  to  several  factors.  First,  the  development  of  the 
use  of  the  relative  phase  features  has  been  based  on  intuitive  judgement  ,  and  not 
proven  estimation  theory.  Therefore  the  use  of  the  phase  features  while  making 
logical  sense,  may  be  far  from  the  optimum  use  of  relative  phase  as  a  classification 
feature. 

An  alternate  method  of  utilizing  relative  phase  information  involves  estima¬ 
tion  of  the  unknown  targets  phase.  Given  a  set  of  coherent  bsckscatter  amplitudes 
one  may  use  the  Cross  Correlation  Algorithm  to  provide  a  maximum  likelihood 
estimate  of  the  unknown  target  phase,  thereby  minimizing  the  conditional  risk  for 
a  uniform  cost  function.  The  phase  estimate  may  then  be  used  to  aid  in  classify¬ 
ing  the  target.  Research  has  shown  that  use  of  the  Cross  Correlation  Algorithm 
will  consistently  improve  classification  performance  over  use  of  amplitude  only 
estimates  [8], 


It  is  critical  to  note  these  results  are  highly  dependent  on  the  noise  model 


used  in  the  study.  Given  a  different  noise  model  where  the  noise  corrupting  the 


amplitude  and  phase  features  is  independent,  different  results  would  he  expected. 
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Parameter  Under  Observation:  p 
Channel  Model:  Relative  Phase  Channel 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
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Parameter  Under  Observation:  p 

Channel  Model:  Relative  Phase  Channel 

Azimuth  angle:  45 

Start  Frequency:  20.0  MHz 

Stop  Frequency:  24.0  Mllz 

Number  of  Frequencies:  10 

IP  Distance  Measure:  j du  f. 

P  :  Key 


_ Non-Coherent  Channel 
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Figure  52:  H'  feature  high  4_MHz  band  weighting  study. 
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Parameter  Under  Observation:  p 

Channel  Model:  Relative  Phase  Channel 

Azimuth  angle:  30 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  5 
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Figure  53:  W  feature  low  <1  MHz  band  weighting  study. 
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Parameter  Under  Observation:  Individual  Distance  Weightings 

Channel  Model:  Relative  Phase  Channel 
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Figure  54:  IF  feature  2  MHz  band  distance  weight  study. 
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Parameter  Under  Observation:  Individual  Distance  Weightings 

Channel  Model:  Relative  Phase  Channel 

Azimuth  angle:  45 

Start  Frequency:  20.0  MHz 

Stop  Frequency:  24.0  MHz 

Number  of  Frequencies:  10 

VU  Distance  Measure:  2<fu  * 

p  :  0.5 
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Figure  55:  W  feature  high  4  MHz  band  dist  ance  weight  study.  ■ 

100  ) 

* 

A 


PERCENT  M6CLASSFICATI0N 

40.0  MO  600  700  80.0  90.0  100.0 


muyyjuuMUUUM 


W^JiraiWJWWVlHWWHIVfW1 


I 


Parameter  Under  Observation:  Individual  Distnnre  Weightings 

Channel  Model:  Relative  Phase  Channel 

Azimuth  angle:  30 

Start  Frequency:  6.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  15 

JU  Distance  Measure:  2^u,Jt 

p  :  0.5 
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Figure  56:  W  feature  low  4  MHz  band  distance  weight  study. 
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CHAPTER  VII 

FADING  CHANNEL  MODEL  STUDY 

7.1  Implementation 

Multipath  propagation  is  a  common  occurence  in  ionospherically  propagated 
signals.  In  the  radar  problem,  multipath  phenomena  generally  results  in  the  re¬ 
ceived  signal  being  the  result  of  the  two  or  more  signals  adding  in  and  out  of  phase 
at  the  receive  antenna.  One  signal,  which  travels  the  expected  path,  may  be  con¬ 
sidered  deterministic  in  nature.  The  other  signals,  which  travel  undesired  paths, 
may  best  be  described  as  a  stochastic  process  since  the  paths  they  take  through 
the  ionosphere  may  be  considered  less  stable.  As  a  result,  the  propagation  path 
may  best  be  described  as  a  fading  channel  [27]. 

A  fading  model  representing  the  estimation  of  features  from  multipath  prop¬ 
agated  signals  as  shown  in  Eq.  (7.1)  was  implemented  in  the  RSSE  program  [28]. 

4F(/,kAf(/,)  =  aAT(/j)f^(/!)+(1_ft2)Lir(/,)Z(/,)f^(i)  i  =  1,JV(7.1) 

where  is  the  deterministic  component,  (1—  a“)^ 

denotes  the  faded  component,  and  a  is  the  weight  of  the  deterministic  component  . 
The  parameter  Aj'(fi)e^f’)  is  the  high  accuracy  complex  backscatter  component 
of  target  T  at  frequency  /t,  Z(i)  is  the  random  amplitude  of  faded  component 
taken  from  a  correlated  or  uncorrelated  random  process,  and  0^(i)  is  the  random 
phase  of  faded  component  taken  from  a  correlated  or  uncorrelated  random  process. 


Construction  of  the  simulated  feature  vector  is  shown  in  Figure  57.  The 
complex  backscattered  data  are  simply  scaled  by  the  chosen  variable  a,  where 
0  <  q  <  1,  to  construct  the  deterministic  component.  Construction  of  the  faded 
components  are  more  complex.  The  magnitude  of  the  high  accuracy  backscattered 
data  is  first  scaled  by  (1  —  o2)2  to  ensure  the  vector  resulting  from  the  sum  of 
the  deterministic  and  faded  components  has  a  total  power  equal  to  the  test  set 
vector  used.  Next,  the  random  samples  Z(i )  are  constructed  by  first  generating 
samples  from  a  white  zero-mean  Gaussian  random  process  with’  a  variance  (cr2)  of 
one.  These  samples  are  then  passed  through  a  filter  with  an  impulse  response  h(i) 
where  i  =0  to  N  —  1.  The  output  process  has  an  autocorrelation  function  given 
by: 

N- 1 

Rz(k)  =  £  h(i)h(k  +  i),  -N  +  l<k<N-l  (7.2) 

i=o 

as  described  in  [9]. 

The  total  power  in  the  process  {  Z(i)  }  is  then  given  by: 

N- 1 

Rz(  0)=  XXO)2  (7.3) 

1=0 

Therefore,  by  selecting  filter  values  h(i)  such  that  the  total  power  of  the  process 
{  Z(i )  }  is  one,  the  resulting  power  of  the  faded  channel  vector  will  be  equivalent 
to  the  power  of  the  test  vector  used.  The  filter  values  used  were  constructed  from 
the  filter  function  by: 

M*)  = - r  *  =  o,...,w-i  (7.4) 

[r£,c-2«‘]5 

where  h(i)  is  the  filter  function,  and  $  is  the  correlation  coefficient  selected  by  user 
the  so  that  Rz{ 0)  equals  1  independent  of  the  value  N  selected.  By  increasing 
$  from  zero  ({  Z(i)  }  highly  correlated)  to  infinity  ({  Z{i)  }  uncorrelated)  the 
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operator  simulates  various  degrees  of  correlation  between  the  faded  components  of 
a  given  measurement  vector. 

The  faded  components  are  then  weighted  by  the  samples  from  the  random 
process  {  Z(i)  }.  Phase  is  assigned  to  each  random  faded  amplitude  by  taking 
samples  from  an  i.i.d.  process  whose  random  variables  are  uniformly  distributed 
over  the  interval  -tt  to  +tt.  The  resulting  complex  deterministic  components  and 
complex  faded  components  are  then  summed,  producing  the  faded  cl.  vnnel  vec¬ 
tor.  Additive  white  Gaussian  noise  is  then  added  coherently  to  the  faded  channel 
vector  as  described  earlier.  The  simulated  vector  represents  the  output  of  a  faded 
channel  with  additive  Gaussian  noise.  By  applying  the  techniques  described  in 
sections  4.3.1  through  4.3.3  to  the  faded  channel  vector,  variations  of  the  faded 
channel  model  such  as  loss  of  phase  information  or  the  addition  of  an  unknown 
multiplicitive  component  may  be  synthesized.  As  a  result,  three  faded  channel 
models  with  characteristics  analogous  to  those  described  in  4.3.1  through  4.3.3 
and  titled  channel  IF,  channel  2F,  and  channel  3F  respectively  were  investigated. 

Three  studies  were  conducted  utilizing  the  faded  channel  model.  First  the  re¬ 
lationship  between  the  deterministic  component  weight  (a)  and  the  classification 
performance  for  the  three  channels  was  investigated.  Next  the  correlation  coeffi¬ 
cient  $  was  varied  while  maintaining  a  constant  deterministic  component  weight. 
Finally,  direct  comparison  between  the  three  faded  channel  models  was  conducted. 

7.2  Deterministic  Weighting  Study 

Figure  58  shows  the  the  degradation  of  classification  of  channel  FI  as  a  is 
decreased.  Note  that  significant  degradation  of  performance  is  only  found  at  0  =  0.5. 
Figure  59  shows  a  significant  degradation  of  classification  performance  at  o  =  0.75 
and  q  =  0.5  for  channel  F2.  Classification  using  channel  F3  demonstrates  an 
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approximate  10  percent  degradation  in  performance  as  o  is  decreased  from  1.0  to 
0.9,  and  from  0.9  to  0.75,  as  indicated  in  Figure  60.  Tims,  these  figures  indicate  that 
the  performance  of  channel  F2  may  be  severely  degraded  by  a  low  deterministic 
weight. 

Figures  61,  62,  and  63  are  similar  to  the  previous  three  figures  with  the  ex¬ 
ception  of  a  smaller  frequency  band  of  2  MHz.  More  significant  degradation  of 
performance  for  channel  FI  is  noted  at  a=0.75  and  a=0.5  as  noted  in  Figure  61, 
which  was  expected  with  the  smaller  frequency  band.  Figure  62  demonstrates  a 
decrease  in  performance  of  25  %  for  a  given  noise  power  as  a  is  varied  from  1.0 
to  0.9.  Likewise,  Figure  63  shows  the  same  large  decrease  in  the  ability  using 
channel  F3  as  the  deterministic  weight  is  decreased. 

Overall,  the  performance  tends  to  degrade  by  a  much  larger  amount  for  the 
same  changes  in  a  at.  the  2  MHz  band.  Therefore  the  deterministic  weight  becomes 
an  increasingly  important  characteristic  of  the  faded  channel  when  the  frequency 
band  decreases. 

A  third  set-  of  Figures  (64,  65,  and  66)  were  constructed  using  one-half  the 
number  of  data  points  over  the  same  8.0  MHz  to  10.0  MHz  band.  Furt  her  decreases 
in  classification  performance  with  decreased  deterministic  weighting  are  expressed 
for  all  three  channels  as  expected. 
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Parameter  Purler  Observaf ion:  Delerministic  Weight 

Channel  Model:  Coherent  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0.5 
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Figure  58:  Weighting  study  for  coherent  4  MHz  band  faded  channel. 
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Parameter  Under  Observation:  Deterministic  Weight 

Channel  Model:  Multiplicitive  Component  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  Mllz 

Stop  Frequency:  12.0  Mllz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0.5 
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Figure  GO:  Weighting  study  for  multiplicitive  4  MHz  band  faded  channel. 
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Parameter  Under  Observation:  Deterministic  Weight 

Channel  Model:  Coherent  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0.5 


3 

S 

V 

s 


a 


n  '  1  ~r  1  1  i - 1  i  '  r< 

oo  5.o  no  e.0  2<xo  ao  mo 


NOSE  POWER  (DB-SM) 


. 1 4  •  .1*  i<i'«U'kU'ft*lV v«m‘  | 


Parameter  Under  Observation:  Deterministic  Weight 

Channel  Model:  Noil-Coherent  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0.5 
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Figure  62:  Weighting  study  for  non-coherent  10  feature  -  2  MHz  band 
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Parameter  Under  Observation:  Deterministic  Weight 

Channel  Model:  Multiplicitive  Component  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight.:  Key 

Correlation  Coefficient:  0.5 


Figure  63:  Weighting  study  for  multiplicitive  10  feature  -  2  MIIz  band  faded 

channel,  112 
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Parameter  Under  Observation:  Deterministic  Weight 

Channel  Model:  Coherent  Faded  Channel 

Azimuth  angle:  45 
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Parameter  Under  Observation:  Deterministic  Weight 

Channel  Model:  Non-Coherent  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MIIz 

Number  of  Frequencies:  5 

Deterministic  Weight:  Key 

Correlation  Coefficient:  0.5 


_ a  =  1.0 

- a  =  0.9 

—  -  —  a  =  0.75 


Figure  65:  Weighting  study  for  non -coherent-  5  feature  -  2  MHz  band  faded 

channel. 
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Parameter  Under  Observation:  Deterministic  Weight. 

Channel  Model:  Mulliplicitive  Component  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  5 

Deterministic  Weight.:  Key 

Correlation  Coefficient:  0.5 
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Figure  66:  Weighting  study  for  multiplicitive  5  feature  -  2  MHz  band  faded 

channel. 
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7.3  Correlation  Study 


Figures  67,  68,  and  69  demonstrate  the  effect  of  the  correlation  of  the  faded 
components  on  the  performance  of  the  three  faded  channels  over  a  2  MHz  frequency 
band.  In  this  band,  channel  FI  does  not  demonstrate  a  decrease  in  classification 
performance  with  increasing  correlation.  Channel  F2  exhibits  a  4%  decrease  in 
performance  when  faded  components  become  highly  correlated  as  shown  in  Fig¬ 
ure  68.  Channel  F3  again  shows  major  improvement  as  a  classification  feature 
when  the  correlation  increases  as  displayed  in  Figure  69. 

Figure  70  demonstrates  the  decrease  in  classification  performance  of  chan¬ 
nel  FI  resulting  from  an  increased  correlation  of  the  faded  components.  In  Fig¬ 
ure  71,  channel  F2  exhibits  a  significant  change  in  classification  performance  at 
$=0.1,  however  the  overall  performance  is  quite  poor  in  all  cases  due  to  the  low 
weighting  (a  =  0.75)  of  the  deterministic  component.  Figure  72  shows  the  increase 
in  classification  performance  of  channel  3  as  the  faded  components  become  more 
correlated.  The  significant  gain  of  3  to  5  dB-sm  in  noise  power  immunity  indicates 
the  correlation  coefficient  may  be  an  important  parameter  in  a  faded  channel  which 
exhibits  an  unknown  multiplicitive  component. 

Figures  67  through  72  indicate  several  relations  between  classification  perfor¬ 
mance  and  correlation  of  faded  components.  Channel  FI,  which  utilizes  intrinsic 
phase  information  tends  to  show  a  slight  decrease  in  performance  with  an  increase 
in  correlation.  This  may  possibly  due  to  the  corrupted  amplitude  of  the  target 
appearing  consistently  smaller  (or  larger)  at  each  frequency  due  to  the  correlation 
of  the  Z(i)  samples.  As  a  result,  the  target  may  appear  to  be  closer  to  a  smaller 
(or  larger)  target  in  terms  of  the  Nearest  Neighbor  Algorithm,  and  therefore  be 
misclassified. 


116 


The  tendency  of  channel  F2  to  increase  performance  as  the  faded  components 
become  uncorrelated  backs  up  the  previously  developed  theory.  However,  since 
channel  F2  removes  all  phase  information  and  passes  only  the  noisy  absolute  am¬ 
plitude  backscatter  components  to  the  classifier,  one  may  expect  to  see  a  greater 
decrease  in  performance  with  increasing  correlation  than  channel  1.  There  is  no 
strong  indication  of  this  behavior,  however. 

Channel  F3  which  adds  the  unknown  multiplicitive  component  which  must 
be  removed  as  described  in  section  4.3.3  demonstrates  a  significant  increase  in 
performance  with  an  increase  in  the  correlation  of  the  faded  components.  As 
reasoned  earlier,  highly  correlated  faded  components  tend  to  make  the  target  look 
larger  (or  smaller)  across  the  band.  As  a  result,  the  average  RCS  of  the  unknown 
target  will  tend  to  be  larger  (or  smaller).  The  process  of  normalizing  the  corrupt 
backscatter  amplitudes  to  remove  the  unknown  multiplicitive  component  also  tends 
to  remove  the  average  contribution  of  the  faded  component.  As  a  result,  the 
normalized  amplitude  feature  is  more  dependent  on  the  deterministic  component 
of  the  faded  channel,  resulting  in  better  classification  performance. 
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Figure  67:  Correlation  study  for  coherent  2  MHz  band  faded  channel. 
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Channel  Model:  Non-Colierent.  Faded  Channel 

Azimuth  angle:  45 
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Figure  68:  Correlation  study  for  non-cohcrent  2  MHz  band  faded  channel. 
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Parameter  Under  Observation:  Correlation 

Channel  Model:  Multiplicitive  Component  Faded  Channel 

Azimuth  angle:  45 

Start  Frequency:  8.0  MIIz 

Stop  Frequency:  10.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  0.9 

Correlation  Coefficient:  Key 
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Figure  69:  Correlation  study  for  multiplicitive  2  MHz  band  faded  channel 
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Parameter  Under  Observation:  Correlation 

Channel  Model:  Coherent  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  0.75 

Correlation  Coefficient:  Key 


Parameter  Under  Observation:  Correlation 
Channel  Model:  Non- Coherent  Failed  Channel 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.75 
Correlation  Coefficient:  Key 
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Parameter  Under  Observation:  Correlation 

Channel  Model:  Multiplicative  Component  Faded  Channel 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  0.75 

Correlation  Coefficient:  Key 
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Figure  72:  Correlation  study  for  multiplicitive  4  MHz  band  faded  channel 
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7.4  Comparison  of  Faded  Channels  Performance 

A  final  study  was  conducted  to  directly  compare  the  performance  of  (lie 
three  channels  for  different  correlation  coefficients,  deterministic  weights,  frequency 
bands,  and  feature  number  combinations.  Figures  73,  74,  and  75  directly  compare 
the  three  feature  types  using  a  2  MHz  band  for  different  values  of  o.  Although  per¬ 
formance  is  relatively  poor  for  both  channels  Fl  and  F2,  channel  F2  significantly 
outperforms  channel  Fl  for  a  —  0.75  and  a  =  0.5  (Figures  74  and  75  respectively). 
Channel  Fl  significantly  outperforms  the  other  two  for  all  values  of  a  as  shown. 

A  comparison  of  the  three  feature  types  for  various  values  of  is  given  in 
Figures  76,  77,  and  78.  Figures  76  and  77  demonstrate  that  for  a  small  ^  (highly 
correlated  faded  components),  channel  F3  may  actually  outperform  channel  F2. 
The  figures  also  demonstrate  the  strong  performance  of  channel  Fl  irregardless  of 
the  correlation  between  faded  components. 

Figures  79  and  80  compare  the  three  faded  channels  for  two  different  values  of 
o.  Both  figures  demonstrate  the  superiority  of  channel  F3  over  channel  F2  in  terms 
of  classification  performance.  In  fact,  the  performance  of  channel  F3  approaches 
that  of  channel  Fl  for  low  noise  power  in  Figure  80. 

Figures  81,  82,  and  83  demonstrate  that  the  resulting  performance  from  using 
channel  F3  features  may  be  significantly  greater  than  using  channel  F2  features 
irregardless  of  the  correlation  of  the  faded  components.  These  figures  also  reinforce 
the  strong  performance  of  channels  utilizing  phase  information.  It  is  important  to 
note,  however,  that  classification  performance  of  a  channel  with  an  unknown  mul- 
tiplicitive  component  (channel  F3)  may  approach  the  performance  of  a  channel 
which  maintains  absolute  amplitude  as  well  as  target  phase  information  (chan¬ 
nel  Fl).  As  the  faded  components  become  increasingly  uncorrelated  in  Figures  81 
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and  83  respectively,  the  performance  of  channel  F3  degrades,  but  still  provides 
relatively  good  classification  performance  at  low  noi&e  powers. 

Direct  comparison  of  the  three  faded  channels  in  term  of  classification  per¬ 
formance  resulted  in  several  interesting  results.  Most  importantly,  the  superior 
performance  of  channel  FI  indicates  that  phase  information  concerning  the  target 
is  maintained  even  for  relatively  small  deterministic  weights  (c*  =  0.5,0.75).  There¬ 
fore,  a  system  which  can  accurately  estimate  intrinsic  target  phase  is  still  desired 
even  if  the  propagation  path  is  subject  to  a  degree  of  multipath  interference.  Sec¬ 
ond,  the  study  indicates  that  under  certain  frequency  band,  deterministic  weight  , 
and  correlation  conditions,  the  removal  of  the  unknown  multiplicitive  component 
from  the  output  of  channel  F3  provides  a  set  of  features  whose  classificat  ion  perfor¬ 
mance  rivals  that  of  channel  FI,  and  is  better  than  the  performance  of  channel  F2. 
Frequency  bands  of  4  MHz,  deterministic  weights  of  0.75,  and  highly  correlated 
faded  components  tend  to  accentuate  the  normalized  features  compared  to  other 
conditions. 

It  is  important  to  note  that  given  a  propagation  path  modeled  by  these  con¬ 
ditions,  which  are  far  from  ideal  in  the  RTI  sense,  classification  may  be  performed 
using  only  relative  RCS  values.  Furthermore,  having  the  ability  to  estimate  ab¬ 
solute  RCS  would  not  provide  an  increase  in  performance.  In  fact,  the  use  of 
the  corrupted  amplitude  estimates  without  normalization  would  seriously  degrade 
classification  performance  (Figure  81).  In  this  case,  even  the  ability  to  estimate 
target  phase  would  only  decrease  the  misclassification  percentage  by  four  or  five 
percent  at  low  noise  power  levels. 

Overall,  the  study  has  indicated  that  if  one  is  given  a  multipath  propagated 
signal,  there  is  little  need  of  a  calibrated  reference  to  estimate  absolute  RCS. 
Rather,  a  stationary  uncalibrated  reference  may  be  utilized  to  estimate  relative 


RCS  of  the  target  to  be  used  for  classification  purposes. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MIIz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.9 
Correlation  Coefficient:  0.5 
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Figure  73:  Highly  deterministic  2  MHz  band  channel  study 


Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.75 
Correlation  Coefficient:  0.5 
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Figure  74:  Partially  faded  2  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.5 
Correlation  Coefficient:  0.5 
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Figure  75:  Highly  faded  2  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  45 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  0.9 

Correlation  Coefficient:  0.001 
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Figure  76:  Highly  correlated  2  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.9 
Correlation  Coefficient:  0.05 
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Figure  77:  Correlated  2  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  45 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  10.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.9 
Correlation  Coefficient:  oc 
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Figure  78:  Uncorrelated  2  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.75 
Correlation  Coefficient:  0.5 
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Figure  79:  Partially  faded  4  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.5 
Correlation  Coefficient:  0.5 
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Figure  80:  Highly  faded  4  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 

Channel  Model:  Key 

Azimuth  angle:  0 

Start  Frequency:  8.0  MHz 

Stop  Frequency:  12.0  MHz 

Number  of  Frequencies:  10 

Deterministic  Weight:  0.75 

Correlation  Coefficient:  0.1 
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Figure  81:  Correlated  4  MHz  band  channel  stud} 


Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MHz 
Number  of  Frequencies:  10 
Deterministic  Weight  :  0.75 
Correlation  Coefficient:  2 


NOISE  POWER  (06-SM) 

_ Coherent  Faded  Channel 

-  -  Non-Coherent.  Faded  Channel 
—  Multiplicitive  Component  Faded  Channel 


Figure  82:  Slighly  correlated  4  MHz  band  channel  study. 
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Parameter  Under  Observation:  Channel 
Channel  Model:  Key 
Azimuth  angle:  0 
Start  Frequency:  8.0  MHz 
Stop  Frequency:  12.0  MIlz 
Number  of  Frequencies:  10 
Deterministic  Weight:  0.75 
Correlation  Coefficient:  oo 
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Figure  83:  Uncorrelated  4  MHz  band  channel  study. 
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CHAPTER  VIII 


CONCLUSIONS 


A  study  of  existing  IIF  radar  systems  with  respect  to  use  as  a  measurement 
facility  has  been  conducted.  Examination  of  a  HF  skywave  radar  system  has  led 
to  the  incorporation  of  two  classes  of  channel  models,  the  Gaussian  and  fading 
channel  models.  Further  variations  of  the  channels  within  each  class  simulated  the 
possible  restrictions  which  the  external  environment  and  propagation  path  might 
place  on  the  measurement  system  in  terms  of  the  estimation  of  various  target 
features.  As  a  result.,  relationship  between  the  various  channels  and  classification 
performance  in  terms  of  various  radar  system  parameters  were  found. 

The  simulation  studies  have  demonstrated  that  classification  performance  is 
largely  dependent  on  the  size  of  the  frequency  band  for  all  channel  models  im¬ 
plemented.  Moreover,  the  results  have  indicated  that  when  the  frequency  band 
is  extremely  limited  (  <  2  MHz),  classification  performance  may  be  improved  by 
increasing  the  number  of  frequencies  at  which  target  features  are  estimated.  Ther- 
fore,  if  one  is  restricted  to  a  small  frequency  band,  performance  may  be  increased 
at  the  cost  of  additional  processing  time. 

The  fading  channel  models  simulated  the  effect  of  multipath  propagation  on 
the  RTI  system.  Results  indicated  that  without  phase  information,  the  estimation 
of  absolute  RCS  feature  may  not  lead  to  additional  classificat  ion  performance  over 
the  estimation  of  relative  RCS  features.  Therfore,  there  would  be  no  advantage  in 
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having  a  calibrated  reference  which  would  allow  the  estimation  of  absolute  IK'S. 

A  variety  of  channels  have  been  simulated.  In  order  to  further  evaluate  classi¬ 
fication  performance,  data  from  existing  HF  radar  systems  is  desired.  Given  a  large 
set  of  Doppler  spectral  estimates  containing  familiar  aircraft,  several  techniques 
may  be  investigated  to  estimate  various  features  of  the  targets.  A  comparison  of 
these  estimates  would  provide  insight  into  the  statistical  properties  of  real  esti¬ 
mates.  In  addition,  a  comparison  of  these  estimates  with  those  simulated  from 
the  EiectroScience  Laboratory  data  base  would  provide  a  measure  of  how  accu¬ 
rately  the  backscatter  coefficients  measured  on  the  compact  range  represent  those 
measured  in  real  HF  radar  system. 
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